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CN Abstract 

The Aharonov-Bohm effect is a fundamental issue in physics. It describes the physically impor- 
tant electromagnetic quantities in quantum mechanics. Its experimental verification constitutes 
a test of the theory of quantum mechanics itself. The remarkable experiments of Tonomura 
et al. ["Observation of Aharonov-Bohm effect by electron holography," Phys. Rev. Lett. 48, 
1443 (1982), "Evidence for Aharonov-Bohm effect with magnetic field completely shielded from 
Qh electron wave", Phys. Rev. Lett. 56, 792 (1986)] arc widely considered as the only experi- 

mental evidence of the physical existence of the Aharonov-Bohm effect. Here we give the first 
rigorous proof that the classical Ansatz of Aharonov and Bohm of 1959 ["Significance of elec- 
tromagnetic potentials in the quantum theory," Phys. Rev. 115, 485 (1959)], that was tested 
by Tonomura et al. , is a good approximation to the exact solution to the Schrodinger equation. 
This also proves that the electron, that is represented by the exact solution, is not accelerated, 
in agreement with the recent experiment of Caprez et al. in 2007 ["Macroscopic test of the 
Aharonov-Bohm effect," Phys. Rev. Lett. 99, 210401 (2007)], that shows that the results of 
the Tonomura et al. experiments can not be explained by the action of a force. Under the 
assumption that the incoming free electron is a gaussian wave packet, we estimate the exact 
solution to the Schrodinger equation for all times. We provide a rigorous, quantitative error 
bound for the difference in norm between the exact solution and the Aharonov-Bohm Ansatz. 
Our bound is uniform in time. We also prove that on the gaussian asymptotic state the scat- 
tering operator is given by a constant phase shift, up to a quantitative error bound that we 
provide. Our results show that for intermediate size electron wave packets, smaller than the 
ones used in the Tonomura et al. experiments, quantum mechanics predicts the results observed 
•*h by Tonomura et al. with an error bound smaller than 10 _ ". It would be quite interesting to 

perform experiments with electron wave packets of intermediate size. Furthermore, we provide 
a physical interpretation of our error bound. 
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1 Introduction 



In classical electrodynamics the force produced by a magnetic field on a charged particle is given by the Lorentz force, 
F = qvx B, where q and v are, respectively, the charge and the velocity of the particle, and B is the magnetic field. In 
regions where the magnetic field is zero the Lorentz force is zero and the particle travels in a straight line. In particular, 
the dynamics of a classical particle is unaffected by magnetic fields enclosed in regions that are not accessible to the 
particle. This also means that in classical electrodynamics the relevant physical quantity is the magnetic field and 
that the magnetic potentials are only a convenient mathematical tool. 

The situation is different in quantum mechanics, where the dynamics is described by the Schrodinger equation that 
can not be formulated directly in terms of the magnetic field. It is required to introduce the magnetic potential. It was 
pointed out by Aharonov and Bohm 21 that this implies that in quantum mechanics the magnetic potentials have a 
real physical significance. Aharonov and Bohm [2.: proposed an experiment to confirm the theoretical prediction. They 
suggested to use a thin, straight solenoid, centered at the origin and with axis in the vertical direction. They supposed 
that the magnetic field was essentially confined to the solenoid. They advised to employ a coherent electron wave 
packet that splits in two parts, each one going trough one side of the solenoid. Both wave packets should be brought 
together behind the solenoid, to create an interference pattern due to the difference in phase in the wave function of 
each part of the wave packet, produced by the magnetic field enclosed inside the solenoid. Actually, the existence 
of this interference pattern was first predicted by Franz [9]. The Aharonov- Bohm effect plays a prominent role in 
fundamental physics, among other reasons, because it describes the physically important electromagnetic quantities 
in quantum mechanics, and since it is a quantum mechanical effect, the verification of its existence constitutes a test 
of the validity of the theory of quantum mechanics itself. 

The case of a solenoid has been extensively studied from the theoretical and experimental points of view. The 
theoretical analysis is reduced to a two dimensional problem after making the assumption that the solenoid is infinite. 
Nevertheless, experimentally it is impossible to have an infinite solenoid and, therefore, the magnetic field can not be 
completely confined into the solenoid. The leakage of the magnetic field was a highly controversial point. To avoid 
this problem it was suggested to use a toroidal magnet, that can contain a magnetic field inside without a leak. The 
experiments with toroidal magnets where carried over by Tonomura et al. |17l 1251 [26] . In remarkable experiments 
they were able to superimpose behind the magnet an electron wave packet that traveled inside the hole of the magnet 
with another electron wave packet that traveled outside the magnet, and they measured the phase shift produced by 
the magnetic flux enclosed in the magnet, giving a strong evidence of the existence of the Aharonov-Bohm effect. In 
fact, the Tonomura et al. experiments [T71 [551 US] are widely considered as the only experimental evidence of the 
existence of the Aharonov-Bohm effect. 

In the case of toroidal magnets, several Ansatze have been provided for the solution to the Schrodinger equation 
and for the scattering matrix without giving error bound estimates for the difference, respectively, between the exact 
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solution and the exact scattering matrix, and the Ansatze. Most of these works are qualitative, although some of 
them give numerical values for their Ansatze. Methods like, Fraunhofer diffraction, first-order Born and high-energy 
approximations, Feynman path integrals and the Kirchhoff method in optics were used to propose the Ansatze. The 
amount of work related to the Aharonov-Bohm effect is very large. For a review of the literature up to 1989 see [15] 
and |18) . In particular, in [18] there is a detailed discussion of the large controversy -involving over three hundred 
papers- concerning the existence of the Aharonov-Bohm effect. For a recent update of this controversy see [23] [27]. 

The paper [4j presents a discussion of a version of the Aharonov-Bohm Ansatz for an infinite solenoid. For recent 
rigorous work in the case of an infinite solenoid see [141 128] where, among other results, it is proven that in the 
high- velocity limit the scattering operator is given by a constant phase shift, as predicted by Franz [9 and Aharonov 
and Bohm [2j. In |16j rigorous mathematical ground is given for the presence of the magnetic potential in the 
Schrodinger operator describing the Aharonov-Bohm effect in the case of a solenoid. In [IT], a semi-classical analysis 
of the Aharonov-Bohm effect in bound-states in two dimensions is given. For a rigorous mathematical analysis of the 
Aharonov-Bohm effect in three dimensions for toroidal magnets -actually in the general case of handle bodies- see [3] , 
where the high- velocity limit of the scattering operator was evaluated in the case where the direction of the velocity 
is kept fixed as its absolute value goes to infinity. A rigorous error bound was given for the difference between the 
scattering operator and its high-velocity limit for incoming asymptotic states that have small interaction with the 
magnet in the high- velocity limit. The error bound goes to zero as the inverse of the velocity. A detailed analysis 
of the Aharonov-Bohm effect in the case of the Tonomura et al. experiments jTTl [25] [26] was given in [3] , as well as 
other results. The results of [3] give a rigorous qualitative proof that quantum mechanics predicts the interference 
patterns observed in the Tonomura et al. experiments [25] [26] [17] with toroidal magnets. The papers [3j[TH[28], as 
well as this paper, use the method introduced in [8] to estimate the high-velocity limit of solutions to Schrodinger 
equations and of the scattering operator. The papers [21], [22], [29], and [30] study the scattering matrix for potentials 
of Aharonov-Bohm type in the whole space. 

In this paper we give the first rigorous proof that the classical Ansatz of Aharonov and Bohm is a good approx- 
imation to the exact solution of the Schrodinger equation. We provide, for the first time, a rigorous quantitative 
mathematical analysis of the Aharonov-Bohm effect with toroidal magnets under the conditions of the experiments 
of Tonomura et al. [TT] [35] EE]- We assume that the incoming free electron is a gaussian wave packet, what from the 
physical point of view is a reasonable assumption. The technical advantage of using a gaussian wave packet for the 
incoming free electrons is that in this case we know very well the dynamics of the free asymptotic gaussian state, and 
we can carry over the estimates of [5] in a precise manner. We provide a rigorous, simple, quantitative, error bound for 
the difference in norm between the exact solution and the approximate solution given by the Aharonov-Bohm Ansatz. 
Our error bound is uniform in time. We also prove that on the gaussian asymptotic state, the scattering operator is 
given by multiplication by e 1 "^* -where q is the charge of the electron, c is the speed of light, h is Planck's constant, 
and $ is the magnetic flux in a transversal section of the magnet- up to a quantitative error bound, that we provide. 
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Actually, the error bound is the same in the cases of the exact solution and the scattering operator. 

Aharonov and Bohm [5] and Tonomura et al. [TTJ [UJ HI] suggested to split the electron wave packet into the 
part that goes through the hole of the magnet and the part that goes outside. Tonomura et al. observed that an 
image was produced behind the magnet that clearly showed that shadow of the magnet and also the hole and the 
exterior of the magnet. They concluded [25j that this indicates that there was not interference between the part 
of the electron wave packet that went trough the hole and the one that either hit the magnet or traveled outside. 
The part of the wave packet that goes outside the magnet can be taken as the reference wave packet. Therefore, we 
only model the part of the electron wave packet that goes through the hole of the magnet. Using the experimental 
data of Tonomura et al. |17[ I25| 126] we provide lower and upper bounds on the variance of the gaussian state in 
order that the electron wave packet actually goes through the hole. We also rigorously prove that the results of the 
Tonomura et al. experiments [17\ I25| 126] , that were predicted by Aharonov and Bohm, actually follow from quantum 
mechanics. Furthermore, our results show that it would be quite interesting to perform experiments for intermediate 
size electron wave packets (smaller than the ones used in the Tonomura et al. experiments, that where much larger 
than the magnet) that satisfies appropriate lower and upper bounds that we provide. One could as well take a larger 
magnet. In this case, the interaction of the electron wave packet with the magnet is negligible -the probability that 



the electron wave packet interacts with the magnet is smaller than 10~ 199 (See Remark 8.12 and Section 9.2)- and, 
moreover, quantum mechanics predicts the results observed by Tonomura et al. with an error bound smaller than 
10 _ ", in norm. 

Our error bound has a physical interpretation. For small variances, it is due to Heisenberg's uncertainty principle. 
If the variance in configuration space is small, the variance in momentum space is big, and then, the component of 
the momentum transversal to the axis of the magnet is large. In consequence, the opening angle of the electron wave 
packet is large, and there is a large interaction with the magnet. If the variance is large, the opening angle is small, 
but as the electron wave packet is big we have again a large interaction with the magnet. 

It has been claimed that the outcome of the Tonomura et al. experiments [T7j [25j [26] can be explained by the 
action of a force acting on the electron that travels through the hole of the magnet. See, for example, [5j [10] and the 
references quoted there. Such a force would accelerate the electron and it would produce a time delay. In a recent 
crucial experiment Caprez et al. [6] found that the time delay is zero, thus experimentally excluding the explanation 
of the results of the Tonomura et al. experiments by the action of a force. In the Aharonov-Bohm Ansatz the electron 
is not accelerated, it propagates following the free evolution, with the wave function multiplied by a phase. Since, as 
mentioned above, we prove that the Aharonov-Bohm Ansatz approximates the exact solution with an error bound 
uniform in time that can be smaller that 10 _ " in norm, we rigorously prove that quantum mechanics predicts that 
no force acts on the electron, in agreement with the experimental results of Caprez et al. 
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1.1 Tonomura et al. Experiments 

The remarkable experiments of Tonomura at al. [171 1251 I26| are widely considered as the only experimental evidence 
of the physical existence of the Aharonov-Bohm effect. Tonomura et. al. constructed small toroidal magnets such 
that the magnetic field is practically zero outside them. In [26j . the magnets are impenetrable and, furthermore, they 
are covered by super conductive layers that forbid the leakage of magnetic field outside the magnets. We denote by 
K := {(xi, X2, x%) e 1R 3 : < f i < (x\ + x 2 ,) 1 / 2 < ^2, 1 333 1 < h} the magnet (fi is the inner radius, the outer radius 
and 2h is the height), and by B(x) the magnetic field. We suppose that B(x) is zero for x outside the magnet. 

An electron wave packet was sent towards the magnet. It was superimposed behind it with a reference electron 
wave packet to produce the interference pattern. The experiments were set up in such a way that the reference electron 
wave packet was not influenced by the magnet, and that the electron wave packet and the reference electron wave 
packet only interfered behind the magnet, were the interference patterns were formed. The observed interference 
patterns provided a strong evidence of the physical existence of the Aharonov-Bohm effect. 

The electron wave packet was much larger than the magnet. It was 3 micrometers in size in the direction of the 
electron propagation and 20 micrometers in size in a plane perpendicular to the propagation direction |24j . It covered 
the magnet completely. Recall that it was observed that an image was produced behind the magnet that clearly 
showed the shadow of the magnet and also the hole and the exterior of the magnet (see [SJJ US]) an d that it was 
pointed out by Tonomura et al. [251 126) , that this indicates that there was no interference between the part of the 
electron wave packet that went through the hole, and the one that either hit the magnet or traveled outside, because 
of the clear image of the shadow of the magnet [UJ [25]. As mentioned before, we will concentrate our analysis on 
the part of the wave packet that goes through the hole, and we will take it as the electron wave packet itself. It is 
either the part of the electron wave packet that goes trough the hole, or a smaller electron wave packet that really 
goes trough the hole. 

1.2 Aharonov-Bohm Ansatz for the Exact Solution 

At the time of emission, i.e., as t —¥ —00, the electron wave packet is far away from the magnet and it does not interact 
with it, therefore, it can be assumed that it follows the free evolution, 



where Hq is the free Hamiltonian. 



d 

ih—<f){x, t) = H (j)(x, (),iel 3 ,iel. (l.i) 



H := — P 2 . (1.2) 
2M V ; 



M is the mass of the electron and P := —iKV is the momentum operator. We represent the emitted electron wave 
packet by the free evolution of a gaussian wave function, (p v , with velocity v, 
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• M 1 _ x 2 

<p v := e^ v - x cp, where <p := ^p^ e"^, (1.3) 

with variance <r smaller than the inner radius of the magnet. We have chosen the variance transverse to the velocity of 
propagation, v, equal to the longitudinal variance in the direction of propagation. In fact, the size of the longitudinal 
variance is not essential for our arguments and we have chosen it equal to the transversal variance only for simplicity. 
Notice that in the momentum representation, e l ^ v ' x is a translation operator by the vector Mv, what implies that 



the wave function (1.3) is centered at the classical momentum Mv in the momentum representation, 



<p v (p) = 0{p - Mv), 

where for any state represented by the wave function (j)(x) in the configuration representation, the momentum repre- 
sentation is given by the Fourier transform, 

■ P 

By the previous analysis, the electron wave packet is represented at the time of emission by the following gaussian 



wave packet that is a solution to the free Schrodinger equation (1.1) 



Vv,o(M) := e~^ H ° ip v {x). (1.4) 
The (exact) electron wave packet, ip v (x,t), satisfies the interacting Schrodinger equation for all times, 



d , - 

ih—ip v (x, t) = Hipv(x, t), x e A := M 3 \ K, t E R, (1.5) 



where 



H:=H(A):=^(P-hAf (1.6) 
is the Hamiltonian and A — jj^A, where c is the speed of light, q is the charge of the electron, h is Plank's constant, 



and A is a magnetic potential with curL4 = B where B is the magnetic field. We define the Hamiltonian (1.6) in 
L 2 (A) with Dirichlet boundary condition at <9A, i.e. ip — for x £ dA. This is the standard boundary condition that 
corresponds to an impenetrable magnet. It implies that the probability that the electron is at the boundary of the 
magnet is zero. Note that the Dirichlet boundary condition is invariant under gauge transformations. In the case of 
the impenetrable magnet the existence of the Aharonov-Bohm effect is more striking, because in this situation there 
is zero interaction of the electron with the magnetic field inside the magnet. Note, however, that once a magnetic 
potential is chosen the particular self-adjoint boundary condition taken at dA does not play an essential role in our 
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calculations. Furthermore, our results hold also for a penetrable magnet where the interacting Schrodinger equation 



(1.5| is defined in all space. Actually, this later case is slightly simpler because we do not need to work with two 
Hilbert spaces, L 2 (K 3 ) for the free evolution, and L 2 (A) for the interacting evolution, what simplifies the proofs. In 



consequence, the electron wave packet is the unique solution, tp v , to the interacting Schrodinger equation (1.5 1 that is 
asymptotic to the free gaussian wave packet, Vv.Oj as i — > — oo, 



ip v (x, t) w ip v ,o{x, i), t — > — oo. (1.7) 

Aharonov and Bohm [5] proposed an approximate solution to the Schrodinger equation over simply connected 
regions (regions with no holes) where the magnetic field is zero, by a change of gauge formula from the zero vector 
potential. Of course, it is not possible to have a gauge transformation from the zero potential everywhere because 
that would imply that the magnetic flux on a transversal section of the magnet would be zero. Hence, the gauge 
transformation has to be discontinuous somewhere. As mentioned in Section |1.1[ in the case of Tonomura et al. 
[TTl E5"l HZtj] experiments the magnet is a cylindrical torus, K. 

We take as the surface of discontinuity of the gauge transformation 

S := {(x 1 ,x 2 ,x 3 ) S E 3 : {x\ + x\fl 2 > r 2 , x 3 = 0} 
and we define the gauge transformation in the domain, V, given by 



V:=A\S. 

Without loss of generality we can suppose that the support of A is contained on the convex hull of K (see Section 7.1 1 



For every x € T>, and a fixed point Xq in T> with vertical component less than —h, we define the gauge transformation 
as follows, 

J x 

where the integral is over a path in T>. Note that for any x € T> with x 3 > the integration contour has to go 
necessarily through the hole of the magnet. 



For any solution to the Schrodinger equation (1.5), <p(x,t), that stays in T>, Aharonov and Bohm [5] propose that 



the solution is given by the following Ansatz, motivated by the change of gauge formula from the zero vector potential, 

(f>AB(x,t) := e ^A,o(x) e -i|ffo e -iAA,a(x) ^Q). (1.8) 

Note that if the initial state at t = is taken as e - iXA - ^ 4>{x, 0) the Aharonov-Bohm Ansatz is the multiplication of 
the free solution by the Dirac magnetic factor e lXA ' ^ [7]. 
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The Aharonov-Bohm Ansatz is expected to be a good approximation to the exact solution if the electron wave 
packet stays in a connected domain, away from the surface S where the gauge transformation is discontinuous. This 
Aharonov-Bohm Ansatz is valid for solutions whose initial data is given at time equal to zero. 



For the incoming electron wave packet that satisfies (1.7) the initial data is given as time tends to — oo and then, 
the Aharonov-Bohm Ansatz has to be modified. To formulate the appropriate Ansatz we define the wave operators, 



W±{A) := W ± := s- lim e'W J e^* 1 * . 

t— f±oo 

where J is the identification operator from L 2 (]R 3 ) into L 2 (A) given by multiplication by the characteristic function of 
A, i.e., J(f>{x) := xa{ x ) 4>{ x ) where, Xa( x ) = l,x & K,xk{ x ) = 0,xeR 3 \A. It is proved in [3] that the strong limits 
exist and that we can replace the operator J by the operator of multiplication by any smooth characteristic cutoff 
function xi x ) £ C°° such that x{ x ) = ®> x 6 K and x( x ) = 1 f° r x m the- complement of a bounded set that contains 
K on its interior. 

The solution to the Schrodinger equation that is asymptotic to the free solution e~' l i Ho (p v as t —> — oo is given by 



i> v := e-**^ W-tp v . (1.9) 

It satisfies, 



lim ||^v - J o|| = 0. 

t — y — OO 



(1.10) 



Using this fact we prove in Section [7] that the Aharonov-Bohm Ansatz for the exact solution to the Schrodinger 



equation (1.5) with initial data as time tends to — oo is given by, 



il>AB,„(x, *) = e^We^V. (1-11) 
what, again, is the multiplication of the free incoming solution by the Dirac magnetic factor e*^- ^ 0. 

It is expected that if the electron wave packet stays in a connected region of space, away from the surface of 
discontinuity <S, the Aharonov-Bohm Ansazt should be a good approximation to the exact solution, i.e., that, 

l/j v W 1pAB,v (1-12) 

The Aharonov-Bohm Ansatz, ipAB.v, is what is observed in the Tonomura et. al. experiments [T71 E51 121)] : as the 
support of the vector potential A is contained in the convex hull of K, for every x whose vertical component is bigger 
than h, Xa,o( x ) is equal to the constant ^;3>, where $ is the flux of the magnetic field over a transverse section of the 
magnet. Then, for x 3 > h, the Aharonov-Bohm Ansatz is given by 
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*PabA x ) = e i &*e- i $ Ho <p y , x 3 > h. 



(1.13) 



This is exactly what it was observed in the Tonomura et al. experiments |17l [ 
The scattering operator is defined as 



S := W*W— 

For large positive times, when the exact electron wave packet is far away from the magnet, and it is localized in the 
region with large positive X3, it can be again approximated with an outgoing solution to the free Schrodinger equation, 



(1.14) 



such that, 



lim \\tp v - J 4>+, v .o\\ = 0. 

t—too 



(1.15) 

The initial data of the incoming and the outgoing solutions to the free Schrodinger equation are related by the 



scattering operator (see Section 3.1) 



<f+,v = S(p v . 

By equations ( |1.10[ ) and ( |1.12|1.16| the Ahar onov-Bohm Ansatz suggests that 



f+.y = S(p v ~ e 



(1.16) 



(1.17) 



i.e., that on the gaussian asymptotic state, (p v , the scattering operator is given by multiplication by e*^*, to a good 
approximation. This also is precisely what was observed in the Tonomura et al. experiments [17, 25, 26J. Furthermore, 
in the Aharonov-Bohm Ansatz ( |1.11| the electron is not accelerated, it propagates along the free evolution, with the 
wave function multiplied by a phase. This implies that in the Aharonov-Bohm Ansatz no force acts on the electron, 
and hence, it is not accelerated. This is precisely what was observed in the Caprez et al. [6] experiments. 

1.3 The Main Results 

As under the free evolution the electron wave packet is concentrated along the classical trajectory, we can expect 
that if the velocity v -that is directed along the positive vertical axis- is large enough, the exact electron wave packet 
will keep away, for all times, from the surface, 5, where the gauge transformation is discontinuous. In consequence, 



the Aharonov-Bohm Ansatz should be a good approximation, and equations (1.121 and (1.17) should hold. In the 



following theorem (see also Theorem 8.10) we prove that this is true under the conditions of the Tonomura et al. 



experiments |17l [213, 26J, provided that appropriate, quantitative, lower and upper bounds on the variance, a, of the 
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gaussian wave function arc satisfied. The requirement for the variance a to lie within the interval below assures that 
interaction of the electron with the magnet and the surface S is small. 

THEOREM 1.1. Aharonov-Bohm Ansatz, Scattering Operator and Tonomura et al. Experiments 
Suppose that the magnets and energies are the ones of the experiments of Tonomura et al.. Then, for every gaus- 
sian wave function, tp, with variance a € [t^Ti ^j*] an d every t € R, the solution to the Schrddinger equation, 
e - l Ti H ( A ) W- (p v , that behaves as e~ l n H ° tp v as t —t — oo is given at the time t by 

<Pab, v := ^^^e-^'V,, (1.18) 



up to the following error, 



\\ e - l i H W-(A)ip v - e lXA -°e- l i Ho ip v \\ < 



7e~^ + 177 x 10 3 e-34^-r^ + 10" 
where, m := M/h. Furthermore, the scattering operator satisfies 



7e~^ + 177 x K^e-i 1 ^^ + lO" 100 . 



(1.19) 



(1.20) 



The main factors that produce the error bound in equation (1.19 1.201 are the terms 



Size of the electron wave packet factor, 



e 2 <* A 



(1.21) 



• Opening angle of the electron wave packet factor, 



33 jamr)- 

e 34 2 



(1.22) 



When the variance a is close to the inner radius of the magnet (the electron wave packet is big), (1.211 is close to 1 



and ( 1.22 1 is extremely small (because in this case o~mv is big ). Then, when the electron wave packet is big compared 



to the inner radius, (1.211 is the important term, what justifies our name. When the variance is small (such that 



amy is close to 1) the factor (1.22 1 is close to one and (1.21) is extremely small ( — is big) and so, the important 



factor is (1.22). Note that when the variance in position, er, is small, by Heisenberg uncertainly principle the variance 



in momentum is big. In particular, the transversal component of momentum is large and the electron wave packet 
spreads a lot as it propagates, what makes the opening angle of the electron wave packet large . This justifies the 



name that we give to ( 1.22 ). Note that in both cases the part of the electron wave packet that hits the obstacle is big. 



When a is big, because the wave packet is big, and when a is small, because the opening angle is big, and even if the 
wave packet was initially small, it spreads rapidly as it propagates inside the magnet and, in consequence, a large part 
of the wave packet hits the obstacle. For variances, a, that are neither to small nor too big the part of the electron 
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wave packet that hits the obstacle is small and the error is very small. In Section [9] we discuss in detail the physical 
interpretation of our error bound and we present a detailed quantitative analysis for a large range of a. 

In particular, we give a rigorous proof that if 1.1592 x 10~ 9 < a < 7.7955 x 10 -6 the error bound is smaller than 
10 _ ". As mentioned above, it would be quite interesting to perform an experiment with electron wave packets that 
satisfy our bounds. One could as well take a larger magnet. In this case the probability that the electron wave packet 



interacts with the magnet is smaller than 10~ 199 (See Remark 8.12 and Section 9.2 1, and quantum mechanics predicts 
with a very small error bound the interference fringes observed in the experiments of Tonomura et al. |171 1251 126) , and 
the absence of a force on the electron, as observed in the Caprez et al. experiment [6 . 

The paper is organized as follows. In Section [2] we introduce notations and definitions that we use along the 
paper. In Section 3 we study the time evolution of the electron wave packet. We define the wave and the scattering 
operators, and we introduce the solutions to the Schrodinger equation with initial condition as time goes to — oo. We 
estimate the solution to the Schrodinger equation when it is incoming, interacting, and outgoing. In Section 4 we use 
the freedom that we have in the selection of the magnetic field, the magnetic potential and the smooth characteristic 
cutoff function to make a choice that is convenient for the computation of the error bounds. In Section 5 we make a 
choice of the free parameters under the experimental conditions of Tonomura et al. |17j . In Section [6] we continue our 
study of the time evolution of the electron wave packet when it is incoming, interacting, and outgoing. In Section [7] 
we consider the Aharonov-Bohm Ansatz for initial data at time zero and for initial data at time — oo. In Section [8] 
we estimate the difference between the exact solution to the Schrodinger equation and the Aharonov-Bohm Ansatz as 
the electron is incoming, interacting, and outgoing. In particular, in Theorem |8.11| we prove our main result that is 
quoted as Theorem |1.1| in the Introduction. In Section [9] give a detailed analysis of the physical interpretation of our 
error bound with quantitative results. In Section [10] we give the conclusions of our paper. In appendix A we prove 
estimates for the free evolution of gaussian states that we use in our work. In Appendix B we prove upper bounds for 
integrals that we need to compute our error bound. 

2 Notations and Definitions 

In this section we collect notations and definitions that are used along the paper. 
The magnet K - see Section [lT]- is defined by the following formula, 

K := {( Xl ,x 2 ,x 3 ) e M 3 : < f i < (x\ + ^) 1/2 < r 2 , |ar 3 | < h) . (2.1) 

We call D the convex hull of K . We use the notation, 

A:=R 3 \K. (2.2) 
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We employ the symbol x — x( x ) = x( x > a ) f° r a twice continuously diffcrcntiable cut-off function that depends on 
the variance of the wave packet, a, - see (1.3). The support of 1 — x is contained in the set 



K := K(a) := {( Xl ,x 2) x 3 ) el 3 :0<n< (xj + xj) 1 / 2 < r 2 , \x 3 \ < h(a)} , (2.3) 

where r\ and r 2 are some positive numbers such that r% < r±, r 2 > f 2 , f\ — r± = r 2 — f 2 and h = h(a) : R + — > K + is 
an increasing function such that h(a) > h for all a in M + . We will write either h or h(a) for the same object. 
We designate by 

e:=fx-rx = r 2 -f 2 , 5(a) := 5 := h(a) - h, (2.4) 
and by D :— D(a) the convex hull of K. 

For every £, u), a € such that < cD < amv, we denote by Z£,, CT (C) the unique solution of the equation, 

(zaAO-C h 4 2 2 77x2u72 = ^~ 1 ' ( 2 - 5 ) 



and for every cti,ct 2 € (0, ri) (see (2.3)) we define 



(ri<Jimv) 2 
ie{i,2} 1 af(mv) 2 + X 2 



Zu>,<T lt (r a (0 ■= Daax(zc 1<71 (0i«£i,ff3(C))) »Vi,<T2 : = mmJA > : 4 \ 2 , , 2 = 1}. (2.6) 



For every a € K + , we define 



,,. /34 a/2000 

^ — V2000) v 33 mv 



mm I a / H amv 



Note that (see equation (11.23) in Appendix A ) 

z(a) > h(a). (2.8) 

For every a £ K+ and every z, s £ K we use the following notation, 

p = p{z) =p{a,z) := , , 2.9 

(a^m^v^ + z^) 1 /^ 



amv 



(a,z,s,Q ■= (C - s) . — jtt-^, 9 inv (a,z) := 9 mv (a, z, z,h(a)), 
(a^m z v 2 + z l yi A 



and 



(2.10) 



(<T,*,Z,C) 2 

T(a,z,s,0~ e- T dr, T(a, z) := T(a, z 7 z, h(a)), (2.11) 

Je inv (a;z,s,-C) 
r e mv (a,z,z,C) 2 

G(a,z,s,() := t 2 e~ T dr, Q(a, z) := Q(a, z, z,h(a)). (2.12) 

Je inv {a,z,s,-Q 

We utilize the symbols fi., c, M and g for the Planck constant, the speed of light and the mass and charge of the 
electron, respectively. We define, 
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M 



We denote by v 6 M' ? the velocity - see (1.3) - and we designate by V := |v|, and v := v/v, respectively, the modulus 
and the direction of the velocity. We suppose that v = (0, 0, 1). We designate by p := — iV x . The momentum operator 
is P := Hp. 

We use the letters B and A for the magnetic field and the magnetic potential, respectively. The details of the 
distribution of the magnetic field inside K are not relevant for the dynamics of the electron that propagates outside 
K, as long as B is contained inside K. Actually, what is relevant is the flux of B along a transversal section of K 
modulo 2ir. See [5] for this issue. We use this freedom to choose B and A in a technically convenient way Then, unless 
we specify something else, we assume that the support of B is contained in K, that the support of A is contained in 
the convex hull of K (what is always possible), and that both are continuously differentiable. In Section [ZJ for any 
given flux in the transversal section of the magnet we explicitly construct a magnetic field and a magnetic potential 
that satisfy our assumptions. We define A :— 4- A, B := 4zB, and 



t]{x,t):= (v x B)(x + pv) dp. (2.13) 
Jo 

We denote by $ the flux of the magnetic field B over a transversal section (TS) of the magnet, 



$ := / B. (2.14) 
Jts 



Then, the flux of B over a transversal section of the magnet is given by, 

$ := / B = — <!. (2.15) 
Jts nc 

By Stokes theorem, for every x — {x\, x%, x$) £ K 3 such that \J x\ + x\ < f\ we have that, 

/oo />oo 
v • A(x + rv)dr, $ = / v • A{x + rv)dr. (2-16) 
-oo J —oo 

Given a function F with domain T> C M. n ,n = 1, 2, • • • that takes values on a normed space C with norm || • ||, we 
denote by Halloo := ess sup{||/(a;)|| : x € V}. 

The vector M = M{\, A, v) = {M 1 (x, A, v), • • • , M 5 (x, A, v)) := (Mi, • • • , M 5 ) G M 5 is given by (v and A and \ 
are defined above in this section ), 
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Mi := ||p 2 xlU + ||XP • AHe., + ||2(p X ) ■ + ||x^ 2 ||cc, 
M 2 := pipxjWoo + PxAU, 
M 3 = \\(px) ■ "lloo + WxA ■ vWoo, 

M 4 := \\x(x)(p ■ A)(x + tfl)||oo + \\x(x)A 2 (x + tv)^ + 2\\A(x + tv) • (px)(«)||oo + 2\\ X (x)A(x + tv)) ■ r](x, t)^, 
M 5 :=2\\x(x)A{x + tv)\\ 0o . 

(2.17) 

The norms in M4 and M5 are taken with respect to x £ R 3 and t € K. 

We define the linear function A : K. —¥ R 5 by the following: given a vector to :— (wi, ■ ■ ■ W5) £ M 5 , we take 
A(w) := (A(w)i, ■■■ , A(w) 5 ) as, 

(2.18) 

The symbols used on the formulae below where defined in this section. Given S\, v £ M + , w £ R 5 and j £ {—00, 0, 00}, 
we define the function Aj, v = A{, : K x R + — > R by 

A- 00 (z, ( t) = A- 00 (z, ( 7) := 

max^,^)^ +max(z,5 1 )- 1 /2(/ lr 2 (CTTOu)3) i/2^)2 + ^M_ z ^li _ 

A° w (z, a) = A°(z, <t) := A.-°°(z, <r) + z.4(w) 4 + ^^H 5 , 

A~(z, a) = A°°(z, ct) := 3A-°°(z, <r) + zA(w) 4 + J^A(w) 5 . 
We will not make explicit the dependence on Si because it will be fixed in our estimates. Actually, Si is a free 
parameter that we introduce to optimize the error bound for the incoming electron wave packet in Theorem |3.1| We 



fix Si in Section 5.2 Note, furthermore, that A. w °°(z,a) is independent of W3 and of W4. We define it as a function 
of w £ R 5 to simplify the statement of our results. 

We define the following quantities, 

C pp (a) = C pp (a,B, X ) := ^^(H A X |U + 2\\r,(x,t) ■ (px)WII) + ^l|PX(*) ' v|U 
C ps (a) = C ps (a,B, X ) := , P ^ W (||XP • tj(i, i)||oo + llxW^, *)!&), 

C sp (ct) = C sp (ct,S,x) := ffV 4mv (IIPX( x )ll°°)' (2.20) 
C„(<7)=tf„(cr,fl,x) := ^i^toCMjIloo), 

_l/2/ 4,2 2 1 >2\l/2 1 //, 5*\2 _(ctttvu)^__ 

ft(C) = ft(C,Z) = ft(C,Z,A) := II A||oo " (J :„1 +C) e"*^ <^W. 
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3 Time Evolution of the Electron Wave Packet 
3.1 Wave and Scattering Operators 



The Hamiltonian operator (1.6) is self-adjoint when it is defined on the domain D(H) := 1 H 2 (A) D Hi.o(A), where 
by H s (A),s = 1,2, •• • we denote the Sobolev spaces and by Hifl(A) we denote the closure in the norm of H\(A) 
of the set (A) of all infinitely differentiable functions with compact support in A pQ . Note that as the functions 
in Hifi(A) vanish in trace sense at dA, H is the positive self-adjoint realization in L 2 (A) of the formal differential 



operator si? (P — HA) 2 with Dirichlet boundary condition at the boundary of A [T31 [TS] - The free Hamiltonian (1.2) 
is self-adjoint when it is defined on the domain D(H ) := ^(R 3 )- Let J be the identification operator from L 2 (R 3 ) 
into L 2 (A) given by multiplication by the characteristic function of A, i.e., 

J(j>(x)= X A(x)4>(x), (3.1) 

where Xa(x) = 1, x G A,\a(x) =0,i£l 3 \A. As mentioned in the introduction, the wave operators are defined as 
follows [20], 

W±(A) = W±:=s- lim e'^Jr'i ff ». (3.2) 



It is proved in [3] that the strong limits (3.2) exist, that they are partially isometric, and that we can replace J by the 
operator of multiplication by any smooth characteristic function, xi x ) G C 2 such that x{x) — 0,x G K and x( x ) — 1 
for x in the complement of a bounded set that contains K on its interior. 

W±(A)=W±=s- lim e^ H xe~^ Ho . (3.3) 

t— >±oo 

It is also known [T3] that the wave operators are asymptotically complete, i.e., that the ranges of W± are the same, 
and that they coincide with the subspace of absolute continuity of H. Moreover, the W± are unitary from L 2 (R 3 ) 
onto the subspace of absolute continuity of H, and they satisfy the intertwining relations, 

e - l i H W ± = W ± e~^ Ha . (3.4) 
Recall that the scattering operator is defined as [2D] , 

S:=WIW^. (3.5) 
3.2 Initial Conditions at Minus Infinity 

In scattering experiments we know the wave packet of the electron at the emission time. Thus, if we want to know the 



evolution of the emitted electron for all times, we have to solve the interacting Schrodinger equation (1.5) with initial 



conditions at minus infinity. As mentioned in the introduction this is accomplished with wave operator W-. The 
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incoming electron wave packet is described at the time of emission (t — > — oo) by a solution to the free Schrbdinger 

(3.6) 



equation, ( 1.1 1, 



As e~ lt yi H is unitary, for all <f>- G L 2 (R 3 ) 



lim 

t— >±oo 



e~ l i H W±<f>- - Je- l i 



Ho 



Then, the solution to (1.5) that behaves as (3.6) as t — > — oo is given by, 



(3.7) 



And, moreover, 



(3.8) 



lim 

t— >oo 



= 0, where 0+ := W*W- (f>- (3.9) 
This means that -as to be expected- for large positive times, when the exact electron wave packet is far away from the 



magnet, it behaves as the outgoing solution to the free Schrbdinger equation (1.1) 



(3.10) 



where the data at t = of the incoming and the outgoing free wave packets (3.6 3.10) are related by the scattering 
operator, 



4>+ = S<f>- 



3.3 The Incoming Electron Wave Packet 



We first introduce concepts that will be used latter in our estimates. 



We define the re-scaled boosted Hamiltonians [3j 128] as follows (see (1.2), (1.6)) 



H 1 = H 1 {w) :-- 



hv 



" v x //„< " ,,v ". H 2 = H 2 (A,v) 



hv 



'■ X H{A) e " nv -' x . 



Recall that m = ^ and v is the velocity (see (1.3)). Let us denote by 



W ± v := e - lmv - x W± e" nv - x 



the boosted wave operators. We have that, 



T^±, v = s lim e lCH2 X {x) e- lCHl 
C— >±oo 



(3.11) 



(3.12) 



(3.13) 
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where ( represents the classical x%— coordinate of the electron at the time t = (/v. 
We notice that, 



-i(H 2 _ -imv-x —i-^H(A) imv-x — iffli _ -imv-x -i-^Ho imv-x _ „-i5^j(P+ mv ) 2 



(3.14) 



The following theorem gives us an estimate of the exact electron wave packet e 1 ^ H W-(A)if v for distances 
Z < — z(cr) < —h(a), i.e., where it is incoming. 

THEOREM 3.1. Let w = (wi, ■••,w&) el 5 be such that un > M L {\,A,w) for i G {1,2,3}. Assume that amv > 
y/M/33. Then, for any Z £E K + such that Z 2i z(&) > h(a), 



(3.15) 



Proof: First we prove (3.151 for W+{A). By Duhamel's formula and (3.14) we have that, 

1 f°° 

\\(W +tV -e tZH2 X z~~ tZHl )v\\<7, / [||mie- l2ffl ^||+2||m 2 -pe-^V||+2mi;||m2-ve- 42ff V||] dz, (3.16) 

zmv J z 

where, 



mi := (p 2 x) - x(P • A) - 2( PX ) • A + A 2 X - 
m 2 ■= (px) - xA. 



(3.17) 
(3.18) 



Equation (3.151 for W+(A) follows from (3.16), Lemmata 11.3 and 11.5 in Appendix A, the facts that the function 
Oinv(o', Z) is decreasing as a function of Z, for Z > 0, that l/w(a) = —9i nv (cr, z(cr)) and the following estimates: 



3/4 



t(z,Si) \a 4 m 2 v 2 + ( 2 

\9 inv (a, Z)\ < amv 



< 2max(z,5 1 )" 1/2 , 



The last inequality follows from the definition of 9 inv (a, Z), since Z ^ z{a) > h(a) (see equation (2.8)) 



We now consider the case of W— (A) . Note that by the uniqueness of the solutions to the Schrodinger equation we 
have that, 



e -iZH 2 (A,v) jjj — e iZH 2 (-A,-v) ^_ 

This is the invariance under time reversal and charge conjugation. Hence, 



W-,- v (-A)iP = W+, V (A)^, 



(3.19) 



(3.20) 



and then, 



(W-- v (-A) - e -iZH 2 (-A,-v) xe .zfl!(-v)j ^ = (w +iV (A) - e iZH *( A > v ) X e-i z HiM) tp . (3.21) 
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It follows that ( |3.15| for W-(-A) and ip_ v follows from ( |3.15| ) for W+(A) and i/j v , and the fact that M(x,A,v) 
M( X ,-A,-v). 



Let L 



be defined as = —a;, for x £ R 3 . Note that, 



(3.22) 



Equation (3.22) implies that, 



(W_ iV (A)y>) o L = W-- v (-A oL)(pL). 
where we used that as = 1 for x in the complement of a bounded set, 



s- lim - = 0. 



(3.23) 



We obtain (3.15) for W-(A) and (p v from (3.15) for VF_(— Ao L), ip- v , x° L instead of \i an d Bo L instead of B 



using equations (3.22 3.23) and observing that M(x,A,v) = M{\ L, —A o L, — v). For this purpose we use B o L 



instead of B in the definition of r\ in (2.13) 



□ 



3.4 The Interacting Electron Wave Packet 



We first introduce an assumption that we use often. 



ASSUMPTION 3.2. Let fa, i € {1,2,3} belong to R+. Suppose that the following conditions hold. 



1. Either fa ^ cr a , i 6 {1, 2, 3}, or fa ^ cr , i € {1, 2, 3}. 



2. Either, fa ^ /i 3 , j G {1, 2}, or fa ^3, i € {1, 2}. 



VFe define /i max := max(/Ji,/j 2 ), /x min := mm(fa,fi 2 ), and iafce f = /Lt min; ^ /i 3 ,i € {1,2} and = ^ max , i/ 

/ii ^/i 3 ,ie {1,2}. FFe denote 6?/ Z := z(/i max ), ^o,ie {1,2,3}; and Z := max l£ { 1:2 }{zi,( Mm „) iMl (/i.(Mmax))}, 
^ a ,i e {1,2,3}. We suppose that Z ^ z y q, Uyfl3 {Ht l mnx)) and rip(fa, z^ v ^ (/i(/x max ))) > 1 for i e {1,2}. 



□ 



The quantities 7 pa , / pp , J aa , and I sp that we use below are defined, respectively, in equations ( 11.26 1, ( 11.32 ), ( 11.35 ), 



and (11.55) in Appendix A. 



LEMMA 3.3. Suppose that Assumption \3.2\ is satisfied and that amv > 1. Then, for every gaussian wave function 
tp with variance a £ [fanin , Mmax] and every £ € R with |£| < z(o~), 
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{ ' ( " ) '° v-A(x+rv)dT) e ~iQH 



< 



(3.24) 



< 



V227! 



-Ml + M; 



(^) 1/2 + ^ 
7r 1 / 4 2o - mt) 



Cpp(o-) 



2 ^PP 



Jpp(Mi 5 M2,M3)+ 



(3.25) 



Cps ^ I pp (^ 1: fl 2j fl 3 ) + Cs ^' T '/ sp (^ 1 ,/i2, /i 3 ) + Cs ^°") / sp (fXi, fj, 2 , fj, 3 ). 

Proof: As in the proof of Lemma 5.6 of j3] (see also [28]) we prove that, 



ELi (fi( x > z ( a ) -t~ z )+ 9i(x, z(a) - C - z) ■ p) 



-izHi 



(3.26) 



where, 



/i0,^) := 2^ [-x(x)(p-^)(x + rv)+ X (a;)(A(a ; + rv)) 2 -2A(a; + Tv)- (p X )(a;) + 
[-X(*)(P ' v)(x,t) +x(x)(r]{x,T))' 2 - (Ax) (a;) - 2tj(x,t) ■ (px)(x)] 



2x(x)A(x + rv) • t)(x,t)] 
1 



2mv 



fa(x) := (px)0) -v, 

9i(x,t) := -—x(x)A(x + tv), 
mv 



9s(x,t) := [-x(x)t](x,t) + (px)(x)} 

mv 



(3.27) 

(3.28) 
(3.29) 
(3.30) 
(3.31) 



It follows that, 

/„* (<r)_C dz \\h{x,z(a) ~ C ~ z)e- tzHl e-< H ^\\ + j'^^ dz\\f 2 (x, z(a) - ( - z) e - lzH i e^" 1 tp\ 



< 



(3.32) 



J z(a) C dz \\h(x)e-^e-^4+J^> z dz \\ 9l (x,z(a)-{-z)-pe 



{<*)-* 



-izH lp -KH 



1 ^ + 



J z(CT) " C dz \\g 2 (x, z{&) -t-z)- pe- lzHl e-^ Hl tp\ 



We estimate the first integral in the right-hand side of (3.32) using equation (11.57), the second using (11.25) and 



(11.31), the third using (11.31), the fourth using (11.59), and the fifth using (11.34) and (11.54). To use (11.59) note 
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that 9 inv (a, z(o~)) = —l/u)(o~). Then, as amv > 1, 9i nv {a, z(er)) 2 > 1/2. After reordering terms we obtain equation 



(3.24). Equation (3.25) is obtained in the same way but using (11.33) instead of (11.31) and (11.56) instead of (11.54) 



□ 



LEMMA 3.4. For Z > ft, 



i J" °° v- A{x+Tv)dr 



e 



X(x)e 



-i/ z -<v-A(x+rv; 



(3.33) 



Proof: By Duhamel's formula and (11.7) 



(x(x)e- i f"*- A ( x+T *'> dT -x{x)e-^o c *-A{*+T*)*r^ e -KB 1(p < \\ x ( x )v ■ A(x + rv) e~ iCBl <p\\ dr 



< 



V2 JZ 



J 7 L 



4„2„2 , ,2 



V2 



(3.34) 



where we used that (ft - r) 2 - (ft - Z) 2 = (r - Z)(r + Z - 2ft). Finally since, (r - Z)(t + Z - 2ft) > (r - Z) 2 , 



< 



Mil oc 



,4„,2„4_i_ r 2 



what proves the lemma. 



(3.35) 



□ 



In the Theorem below we estimate the exact electron wave packet e l "H H W±(A)(p v for distances £ such that, 



ICI < z (°")- As z(cr) > ft(cr) -see equation (2.8)- this is the interaction region 



THEOREM 3.5. Suppose that Assumption \3.S\ is satisfied, and, furthermore that amv > 1. Let w — (u>i, • • • , w§) £ 
K 6e such that W{ > M,(x, A, v) /or i £ {1, • • • ,5}. Then, for every gaussian wave function tp with variance a £ 
[/"mm, A*max] a™d ewer?/ ( £l with \(\ < z(o"), 



\\e^& H W±(A)^-xe^tf~*- A t*^* r e- i & a °<Pv\\<e WA° W (z(a),a) + C pp (a)7 pp ( Ml)M2 , M3 )+ 

^2 -WW.) l"2,M3,±C) +C S p(c r )^p(Ml I A 1 2, ) "3) + %^J as (/il,M2)M3,±C) + |^(±C) 



(3.36) 



||W±(A)^ V - xe-^o^v-Af.+rv^r^n < e " =W (A"°° (2(a), a) + 4fU( W ) 4 + #^H 5 )+ 



C 2^ sp(Ml)M2,M3) + 5^(0, z(cr)). 



Proof: 



(3.37) 



We prove (3.36) for W + (A), the proof for W_(A) follows as in ( 3.19|3.23 ). Note that by the intertwining relations 



of the wave operators (3.4) and by (3.14) we have that, 
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(3.38) 



We use the intertwining relations and (3.14) again to obtain, 



||^ + , v (A)e-< H V-Xe~ l J'o ±oc ^ A ( 2:+T ^ <iT e-^ ff Vll < ||(W + , V (A) - e^M-OHi^iz^-QH^^^^ 

\U e i{z{a)-C)H 2 -i( z ( ff )-C)Hi _ xe -»/ [ f (<T) " C v-A(x+Tv)dT^ e - l C«i ( p|| + 



(3.39) 



Equation (3.36) is obtained by (3.38), (3.39), Theorem 3.1 equation (3.24) and Lemma 3.4 The proof of (3.37) is 



similar, but instead of (3.24) we use (3.25) 



3.5 Estimates for the Scattering Operator 



We first prove the following lemma. 



LEMMA 3.6. Suppose that the conditions of Theorem \ 3. 5\ are satisfied. Then, 



CpP 2 a) Ippfau M2, Ma) + 2 j pp(Mi»M2»M3) + 2 ^p(Mi»M2,M3) + CI " 2 (tT ^ S p(Mi i M2,M3)+ 

|W(0,Z(<7)). 

Proo/: As W+ = /, 



(3.40) 



< 



^ e -i/ <>o v-^+rv)dT x _ W+.v) (y5 + ||(1 - XMI + (e" i( ^ °° v-A(x-frv) dr+*) _ g-i / °° v-A(aH-rv) dr)^ 



(3.41) 



Since J_ v • A(x + rv) c?r = <f> for a; in the cylinder {x ^ M? : x\ + x\ < rf}, (3.40) follows from Theorem 
the following estimates, 



3.5 



and 



\\(l- X (x)M<e- r ^ , 



(e 



-i(J 00 v-A(x+tv) dr+$) _ g — i /" v-A(a;+rv) dr 



< 2e- r i/ 2CT2 



(3.42) 
□ 



In the theorem below we approximate the scattering operator by its high- velocity limit (see [3]). 
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THEOREM 3.7. Suppose that Assumption 3.2 is satisfied. Let w = (wi, • • • , W5) G K 5 be such that Wi > Mj(x, A, v) 
/or i G {1, • • • , 5}. Then, for every gaussian wave function (p with variance a e [^ m i n , ^max]; 



|| (5 - e l * X ) ¥>v|| < 3e"0 + e (2A-°°(z(<7), a) + z(a)^( W ) 4 + $&A(w) 6 ) + C„(o)I„(pn, M2, Ms)- 

C ps (o-)7p P (^i,^2,A i 3) + C sp (o-)/ S j,(^i,^2,M3) + C , ss (CT)/ 5 .p(/ii,/i 2 ,/i3) + 72(0, z(cr)). 



(3.43) 



Proof: We denote, 



We have that, 



g g-imv-i S e' lmv '' 



(3.44) 



||(5-e**x) Vv|| = ||(Sv-e**xH| = |Wf,v (W-.v-xW^'^"^^'" 7 ) ¥> + 



(3.45) 



Equation (3.431 follows from Theorem 3.5 Lemma 3.6 and (3.451 



3.6 The Outgoing Electron Wave Packet 



In the following theorem we estimate the exact electron wave packet e 1 ^ H W-(A)ip v for distances £ in the outgoing 
region, ( > z(a) > h(a). 



THEOREM 3.8. Suppose that Assumption 3.2 is satisfied. Let w — (wi, • • • , w$) £ M 5 be such that Wi > Mi(\, A, v) 
for i e {1, • • • , 5}. Then, for every gaussian wave function ip with variance a G [^ m in, /^max] and every ( el wzi/i 



||e- l A^W_(A)^ v - X e l *e- 8 * H Vv|| <3e-^+e A~(z(<r), a) + C„{a)I„(jn, /i 2 , // 3 )+ 

C ps {<r)Ipp(ni, V2,V3) + C S p(o-)/ S p(/ii,/i 2 ,M3) + C ss (a)I sp (fi 1 ,fi 2 ,^3) + U(0,z(a)). 



Proof: 



Using the definition of S 1 (see (3.5)) and the fact that e is unitary we get, 



(3.46) 



\\ e - i ^ H W-(A)ip v ~ xe^e-^^vll = \\W-(A)ip v - e l ^ H X e i ' i 'e- l ^ Ha ip v \\ < 
He 4 * (w+(A)p v - e'a H xe" , ^ iI °) Vv|| + ||W_(A)</? V - W+(A)e^ ip v \\. 



Furthermore, 



(3.47) 



\\W_(A) Vv -W + (A)e^^\\ < \\W-(A)<p v -W + (A)S<p v \\ + \\W + (A)(S-e**)<p v \ 



(3.48) 
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Since the wave operators are asymptotically complete [13], the operators W± W± are the orthogonal projector onto 
the common range of W±. Then, W + W— — W-, and we have that, 

W-(A)ip v - W+(A)Sip v = W-(A)ip v - W+(A) WUA)W-(A)cp v = 0, 



and by (3.47 3.481 



\e-^ H W-{A)p v - xe l *e-^ H Vv|| < \\W+(A)^ V - e^ H X e~ l ^ H °^ v \\ + \\S<p v - e l *(^ v || 



(3.49) 



The inequality (3.461 follows from Theorems 3.1 and 3.7 and from equation (3.49). 



4 The Magnetic Field, the Magnetic Potential and the Cutoff Function 



We have proven in Theorem 4.1 of [3] that the Hamiltonias (1.6 1 with Dirichlet boundary condition on dA that 
correspond to two different magnetic fields contained inside the magnet, and that have the same flux $ modulo 2tt are 
unitarily equivalent. We have also proven in [3] that the scattering operator only depends on the total flux $ enclosed 
inside the magnet, modulo 2n. This implies that without losing generality we can assume that 



|$| < 2tt, (4.1) 

what we do from now on. This also means that we have a large freedom to choose the magnetic field, as long as it is 
contained inside the magnet. As mentioned in the introduction, we also have a large freedom to choose the smooth 
cutoff function x- We use this freedom to choose the magnetic field, the magnetic potential and the smooth cutoff 
function that is convenient for the computation of the error bounds. Below we construct a magnetic field inspired in 
the experimental results of Tonomura et. al. [25j . We also choose a magnetic potential and a cutoff function, and wc 
provide bounds for them. 

4.1 Mollifiers 

We denote for z e K, 

! f |*| <1, 

m := - \ (4.2) 
0, |z|>l, 



where, 



For e > we define, 



1 e-^^-^dz. (4.3) 



A{z) := U(z/e), (4.4) 

£ 
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and for every a, b G K, with a < b and every e G K + with e < \ (b — a) , we take, 

1pa,b,e( z ) : = 

Then, 



i 1, z G [a + e, b — s], 

dyip £ (z -y) = 

0, z [a — e, b + e]. 



IV'o, 



b.e 



(4.5) 



(4.6) 



b,e\ 



< 



ie e 



b.e 



2N 
- — > 



where JV := 2 e -( 3 / 2 +V 3 / 4 ) (3/2 + V3?4) 2 (l - (3/2 + ^/i/l)- 1 ) 11 ' 2 



(4.7) 
(4.8) 



4.2 The Magnetic Field 



Recall that the magnet is the set, 



K 



:= {(si,a 2> ar 3 ) G M. 3 : < f i < (a? + ^) 1 /2 < ^ | X3 | < fcj 



(4.9) 



We use cylindrical coordinates: for [x\,X2, £3) G M 3 , we take r := (x\ + x^) 1 / 2 , < 8 < 2n, X3. For e < T2 4 Tl , 5 < 
|, we define, 



B = B(x,e,5) := ^fi+efy-eA r ) ^-h+S h-s s( x s)\- sin ^^ cos& ^ )^ 



where for a transverse section of if,TS, 



(4.10) 



C g 5 := / ipf 1+s ,f 2 -e,i(r) ip_~ h+ s h-s s( x s) > 2 ~ 2 <5) - ri - 4e). 



(4.11) 



Then, V ■ B = and the flux of B over any transverse section of K is <E>. 

This choice of -B, that is approximately constant along any transverse section of K and is directed along the unit 
vector (— sm(0), cos(0), 0) is inspired by the experimental results of Tonomura et al. [25]: in Figure 4 (a) of [2"5] . 
the fringes on the shadow of the magnet suggest that the component of the magnetic field that is orthogonal to a 
transverse section of the magnet is constant over this transverse section. 



By (4.6 4.7 4.111, 



dx 



B 



(h - 25) (fa -fx- 4e) 

7T 



,1 1 x 

-(— : + ^), J = 1,2, 



(h - 25) (f 2 - f 1 - 4e) iee r\ 



(4.12) 
(4.13) 
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dx 3 ' 



B 



< 



(Ji-25) (f 2 -fi-4e) LeS' 



(4.14) 



With this choice of B we have that (see ( 2.13 ) ) . 



W^-OIU < 2h- l — -~— ; — . 

(ft - 28) (r 2 - 7"i - 4e) 



(4.15) 



Hp • »?(a;,i")||oo < 2ft -~ 



(ft - 28) (f 2 - fx - 4e) \ iee r i 



1 1 



(4.16) 



4.3 The Magnetic Potential 



The potential A = A(x, e, 8) associated to the field B = B(x, e, 8) satisfies the differential equation V x A = B. As B 
has no vertical component, we can take A parallel to the vertical axis. 

'(vim) j 
where (2/1,2/2) is any point with | (2/1,2/2)! > fi and the line integral is over any curve in K 2 that connects the point 
(2/i;2/2) with (xx,x 2 )- The value of A is independent of the curve chosen. The potential A has support in the convex 



A = A(x,S,S) := ■^;^_i l _; > j,_; i jir: i ) ((>.(). / :.;(/')( cos t>. suit 



hull of K, that we denoted by D. Moreover, by (4.6 4.7) 



l^lloo < - 



dx 



-A 



(ft - 28) (r a - f 1 - 4e) 



< — 



— (^2 - fi), 



(ft - 25) (f 2 - fi - 4e) 



,.7 = 1,2, 



(4.18) 
(4.19) 



dx 3 



.4 



< — 



(ft - 25) (f 2 - fi - 4e) 6e5 



7 (ra - ri). 



(4.20) 



4.4 The Cutoff Function 



We use the freedom that we have in the choice of the cutoff function x( x ) to select it in a convenient way. Take 



< e < fx, 8 > 0. We define (see (2.4)) 



ri := fx — e > 0, r 2 ■= r 2 + e, ft := ft + 5. 



(4.21) 



We define 



X(x) := 1- ?/Vi+ e /2,r 2 -e/2, £ /2(r) 1p-h+S/2,h-5/2,S/2{x 3 ). 



(4.22) 
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Then (see p3|) 



Moreover, by (4.6 4.7 4.8), 



0, x G K, 



l, iel 3 \x. 



Ilxlloo = 



dxi 



< —,3 = 1,2, 

Lee 



< 



Led ' 



We denote by 



P X < 



8iV 



8iV 



lb 2 er\ie lS 2 



I:= I(ft-2«)(f 2 -r 1 -4£) J 



7 r 2 -ri 



We designate by m(x) =m:= (mi(x), ■ ■ • , Ws(x)) G ^ 5 the vector with the following components, 



mi(x) = mx := 



^ + + ^ + ( 2 + (f 2 - f x )~ ) J- 1 + -f J + J 2 



m 2 ( X ) = ra 2 := 2 + J^) + 2J, 
2 

"^(x) = m 3 := — + J, 



1 



m 4 ( X ) = m 4 := ( 2 + (f 2 - f x )-J- ) I^ 1 + J 2 + ±J, 



m s(x) = m 5 ■= 2 J. 
Now we define the following quantities, 



Cpp(er) 

C S p(c) 
c ss (cr) 



8N 2 87V Ah 4 



2/i 



7T ' cjtti^ \ Lee 



1 4/t 



tee r*i 
4 

iSe 



*" 4 ™ \^e 2 tene t6 2 I Lee I n^iSe' 
1 1 



, \ m 5 (cr 4 m 2 u 2 + C 2 ) 1/2 1/9 1 ih-z)*(<, m «f 
R(C,Z) = R(C) := — — 7r 1/2 e 2 ^m^+c* . 
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REMARK 4.1. For the field, the potential and cutoff function constructed in this section we haue that, 
Mi <m u ie {l,---,5}, 

C pp {cr) < Cpp(a), C ps (a) < c ps (a), C ss {a) < c ss (a), C sp (a) < c sp (a), (4-31) 
ll((,Z)<R((,Z). 
Proof: the Remark follows from explicit computation. 

□ 

We introduce some notation that we use below. We define the vectors A J (v, fh) = A J := (A\, A^ 2 , A^, A^^, A^-J. 
for j € {— oo, 0, oo}: 

A-°°(v,fh) = A-°° := (mvn^ + mvr 2 (^ (1 _ 5xl ^ 10)1/a ■ 
134.99ft^f^,-134.99/i^f^,0), 

A~°(v, fh) = A~° := (Ar°°, A v ^, Aq 00 + 135.91/U(m) 4 , Al^ 2 + (4.32) 
135.91/U(m) 5) ^^(1 + 1.11 x l - 6 )V 2 iff /i), 

A°°(v,m) = A°° := (3A^°°, 3A-™, 3Aq 00 + 138/U(m) 4 , 3Al^ 2 + 138/U(m) 5 , 0). 
Finally, for j G {— oo, 0, oo} we denote, 

A J (cr,w,m) = A 3 {a) := ^ A^V. (4.33) 

iS{l, 1/2, 0,-1/2,-1} 

5 Tonomura et al. Experiments. Continued 
5.1 Experimental Data 

We consider the 2 different magnets with their dimensions given in table I of [IT]- We denote them by {^j}jg{i,2}, 

Kj := {x = (x 1 ,x 2 ,x 3 ) e M 3 : h,j < \Jx\+x\ < f 2 ,j, \x 3 \ < h}. (5.1) 

We use the notation 

Xi,j€ {1,2} (5.2) 



for the corresponding cutoff function constructed in Section 4.4 
The height h is 10 _6 cm for both magnets and 

fn = 1.5 x 10~ 4 cm, 
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f 2j i = 2.5 x 10 4 cm, 



h,2 = 1-75 x 10~ 4 cto, 
f 2 ,2 = 2.75 x l(T 4 c77i. 

In the Tonomura et al. experiments |26j the electron has an energy of 150 keV. In this experiments they consider 
impenetrable magnets as we do in this paper. In the experiments |25j they consider penetrable magnets and energies 
of 80 keV, 100 keV and 125 keV. Since our method applies also in the case of penetrable magnets, we will consider 
in our estimates below the two extreme energies and an intermediate energy, although the most important one is the 
one of 150 keV that is the one used for the case of impenetrable magnets. Thus we consider the following energies. 

E l = 150fceV; 
E 2 = 100 keV, 
E 3 = 80 keV. 

They used an electron wave packet that might be represented at the time of emission ( t — > — oo) by the gaussian 
wave function, 

(±) , ' , (4.y' t ->w™ t -w.-£i. (5.3) 

\a 2 z irj \a 2 r TiJ 

The transverse variance of the wave function a r is several times the radius of the torus (r 2j , j = 1,2), so the 
electron wave packet covers the magnet. 

The part of the wave packet that goes through the hole of the torus has a different behavior than the one that goes 
outside the hole. There appears to be no interference between those two parts of the wave packet, because a clear 
figure of the shadow of magnet is formed behind the torus. This was pointed out by Tonomura et al. [35], [55]. We 
can, therefore, model only the part of the electron wave packet that goes trough the hole of the magnet. Hence, we 
take the transverse variance a r smaller than the inner radius of the magnet. The anisotropy of the variance (a z ^ a r ) 
does not introduce new ideas to the analysis and all the proofs that we do assuming that a z = a r can be done in the 
same way if a z ^ a r . We obtain similar results in both situations. Taking a z ^ a r complicates the notations and, 
therefore, for simplicity, we will assume that a z = a r = a. So, when emitted, the electron that goes trough the hole 
is represented by, 

lM*,t) := r^e-^^^-e"^, (5.4) 
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with the variance a smaller than the inner radius of the magnet. 

The real electron wave packet, under the experiment conditions, that behaves as (5.4) when the time goes to — oo 
is given by the wave function (see (1.9)), 



tp v (x,t) := e-^ H W-<p v = e~ i ^ H W^(p v . (5.5) 

Remember that we take v = (0, 0, v) and that £ := vt is the classical position of the electron, in the vertical direction, 
at time t. 

The energy for the free wave packet (or of the perturbed wave packet at — oo) is given by 



When a is big ( amv » 1 ) the second factor is much smaller than the first. If we take for example amv > yl5 
the second factor is less that 1/10 times the first. Therefore, when amv >> 1, we can suppose that the energy is 
given by the classical energy, m v 2 . With this assumption we can calculate the velocities, and the velocities times m 
corresponding to the energies Ei,E 2 ,E 3 : 

Vt = 2.2971 x 10 w cm/s, mv x = 1.9842 x 10 10 cm~\ 

v 2 = 1.8755 x 10 10 cm/s, mv 2 = 1.6201 x 10 10 cm~\ 
v 3 = 1.6775 x 10 10 cm/s, mv 3 = 1.4491 x 10 10 cm -1 . 
For now on we suppose that the obstacle K is either K x or K 2 and that the velocity v is either v\, v 2 or v 3 . 

5.2 Selection of the Parameters 

We have obtained rigorous upper bounds for the difference between the exact solution to the Schrodingcr equation 
and the Aharonov-Bohm Ansatz, and for the difference between the scattering operator and its high-velocity limit. 
These bounds hold for any choice of the parameters Si, S, e, 8 and e. We use this freedom to choose these parameters 
in a convenient way. From now on, we choose the parameter S% > such that 

np(Si) = 1. (5.7) 

This choice is made to optimize the error bound in Theorem |3.1| This theorem was proven using Lemmata |11.3| pi 1 . 5 1 



For example, for the convergence of the integral on the left-hand side of equation (11.10) we need the decay of p{a, z) 
for large z, but for z small this factor is very large. For this reason we split this integral in two regions (where we use 
different estimates) introducing the parameter S±. 
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Furthermore. 



2ll(l 



S-= -i- 

° • 100' 

5 := max(10 cr, h), 



(5.8) 



50 ' 



This selection was obtained using numerical estimates to optimize the error bound for the time evolution of the electron 
wave packet. 

6 The Time Evolution of the Electron Wave Packet. Continued 

LEMMA 6.1. For the data used in the Tonomura et al. experiments, v G {v\,V2,v^\ an d K € {K^^K^, suppose 
that a e \— ,7=1/21 and (eK. Then, 



e 5; ^A^ 0O (z( ( t), ( 7) < e-»^-3 J -A- oo ((r) + 10- 



-420 



e 5oc^A!kCs(tr),ff) + iii(C,z(<T)) < e-H^a^ A°(tr) + 10 



-420 



e A°°(z(a), a) + R(0, z(o-)) < e" if A°° (cr) + 10 



-420 



Proof: 



(6.1) 



First case, cr € [<To, y\- 

As w(cr)" 1 < w||crm«, we have that 



1 < 



(crmt>)^ 



(amv) 2 — w(cr) 2 



< 34. 



For these values of cr, w(cr)- 1 = V2000. Then, using (|ll.23[) and the experimental values we get, 



(6.2) 



2.1023 x 10~ 6 < z(a) < .0673. 



We also have, 



(6.3) 



.0042 < S! < 303.8306. 



Using (6.3) and (6.4) we get, 



(6.4) 



(ftr 2 cr 3 TOV) 1/2 (max(z(cr),S' 1 )) _1/2 < 2.9127 x 10 5 



(6.5) 
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and 



34 x 2000 ' 



(6.6) 



Now we note that (see the definition of i?(C, z{a)) in (4.30)) 



R((,z(a)) < ^ (gW + z(g)2)1/2 ^ e -^ ; mz(a)) < ^/ Ve 



m 5 1/2 _ i^ = £W^ 

TT ' (TP 2ct 2 

2 



(6.7) 



We bound the quantities A J , je {— oo,0, oo} uniformly for a G [ctq, ^f] and for the experimental energies and 



magnets, using (6.3 



6.4 



6.5 



6.7) and the smaller experimental values of f\, (f-z — fi), ft and mv to determine the 

_ i 

components of in. We use the fact that for the values of sigma that we consider, e ^w) 2 < e - 1000 to obtain, 



e 2e( .)2 A-°°(;j:(o-),a-) < 10' 



-420 



e"^A°(z(a),a) + iR(£z(tr)) < lO" 420 , 



(6.8) 



e 5 ^ ? A°°(z(ct), cr) + i?(0,z(cr)) < 10 



-420 



Second case, a G [-r|, ctq]. For these values of cr, 



(crmu) 
(frm?;) 2 +cj(er)~ 



= 34, then by (11.23), 34 ft + V34 V33 ft < z(a) < 



34 ft + V34 J %u A m 2 v 2 + 33 ft 2 and by triangle inequality z{a) < 34ft, + \/33 cr 2 mw + 34ft and then, we have that, 



134.99 ft < z(a) < 136.82 ft. 



(6.9) 



It can be verified that, 



max(z(cr), Si) = Si < amvri, 

max(z( C r),^ 1 )- 1 /2 ( ^2 (J 3 TO 3, ; 3 ) l/2 < (^^(A)^ l_ 10)1/a , 



{a l m 2 v 2 + Z {a) 2 ) 1 ' 2 



< (1.11 X 10~ 6 + 1) 



1/2 136.82/t 



where in the last inequality we used (6.9). Using (6.9) again we get 



(6.10) 



n , . Ill* I/O 1 / 133.99/i \2 ln H 

R(0,z{a)) < ^7r 1/2 cre-5 ( ^^ ) < 10~ 10 



(6.11) 



Finally we obtain (6.1| using (2.19), (4.32), (6.7), (6.9), (6.10), (6.11) and the fact that A±(m) < Ai{m) and 



A^(m) < As{m) (note that in this case e 2i 
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e 34 2 J. 



33 (trmv) 



REMARK 6.2. For j G {— oo,0,oo}, e 34 2 h?{a) is decreasing on the interval [^,00) 



Proof: Calculating the numbers A\ we find that A^ > for i € {1,1/2,-1}, and also A ^^ > 0. The other 
components of the vectors A- 7 are negative. We suppose that j G {— 00,00}, the case j — can be done in the 
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same way (the term A^ 1 , 2 is manipulated as the term A° 1 ). Since A?(<j) > and amv > 4.5, we have that 
£ e -^^^ A j < e -l i ^(-6 1 (cr)+6 2 (CT)), where b x {a) = |f 4.5mw(A{cr+A J 1/2 cr 1 /2) > o and o 2 = -§§ 4.5mv(A J + 
-^■-i/2 cr_1 ^ 2 ) Sie{i 1/2 o -1/2} *A^ct 1 _1 > 0. As 6 X is increasing and 6 2 decreasing, — bi + b 2 is decreasing, as 
-61O + 62O < 0, we have that £ e "^A' < for a e [££,<r ]. 

□ 

Below we introduce a partition of an interval that is adapted to the order of magnitude. 

DEFINITION 6.3. For any number a > we designate by O a G Z the order of a, (i.e. O a is such that 10°" < 
a < 10 Oa+1 ). For an interval [a, b], a > and a positive number No we define the partition V(a,b,N ) := {pi}^ = i 
(p« < Pi+i^i € {1, • • • , fc — 1}) as follows: 

• case 1: b < 10°° +1 . If b - a < A^IO " we take k = 2, Pl = a, p 2 = b. If b - a > N o 10°* we take k > 3, 
Pi = a, pfe = 6 and Pi,i £ {2, • • • k — 1} such that < Pi+i, Pi+i — pi = A^IO " for i e {1, • • • , k — 2} and 
Pk-Pk-i <N 10°«. 

• case 2: b > 10 Oa+1 . For every j £ {0, • • • , Ob — O a } we define a set as follows. We take V° as in the case 1 
but taking 10 Oa+1 instead of b. r p°i>- ^ \ s taken as in the case 1 taking 10 Ob instead of a. If Of, — O a > 2, for 
j G {1, • • • , Ob — O a — 1} we define as in the case 1 taking lO "^ instead of a and 10 Oa+J+1 instead of b. 
Now we define V(a,b,N ) = Uj G { ,- ^^oj^- 

DEFINITION 6.4. We denoie 6y the following sets: 

S i : = niog^^iogfe-OOOS), E 2 := niog^^iogcSyiso-OOOS), E 3 := P( i , g (^, HT 5 , .0008), £ 4 := 
-P(10- 5 ,l.l x 10- 5 ,.0001), £ 5 := x 10" 5 ,1.3 x lO" 5 , .0002), E 6 := V(1.3 x 10" 5 , 1.7 x 10- 5 ,.0004), E 7 := 

P(1.7 x HT 5 , 2 x 10- 5 , .0008), £ 8 := P(2 x 10" 5 , ^, .0015), £ 9 := 7>(10- 6 , log(1 r ^ )250 , 1000), S 10 := V(a , lO" 6 , 1000), 
Sn :=V(^ v ,a ,.l). 

LEMMA 6.5. Suppose that the energies and magnets are the ones used on Tonomura et al. experiments. Let 
Hi e R + , i e {1, 2, 3}. Suppose that {/Zj} 2 =1 is contained in one of the sets Ej for j e {1, • • • , 11}. We take ^3 = 10~ 6 
if {Mi}?=i * s contained in Ej /or j G {1, • • • , 10} and we take ^3 = a n if {/ij} 2 =1 is contained in the last set. We 
suppose furthermore, that [i\ and [i 2 are consecutive numbers in the set where they belong and /ii < /i 2 . Then, for 
every a G \pi, /i 2 ] and every £ S R wii/i |£| < z(ct) we /wroe i/iai, 

Cpp{<?)Ipp(p-l, M2, M3) + ^^4^ ^ G"i '^2, A*3, C) + c sp (cr)/ sp (^1,^2,^3) + ££ 4 £ ^ / s S (Ail,Ai2,M3,C) < 
4e-^ + 10- 3 e-I i ^A ( ( 7) + lO" 101 , 

(6.12) 

Cpp(o-)/pp(Ati, jti 2 ,A*3) + c P s(o-)Ipp(vi, ^2,^3) + c sp {a)I sp {in,ii2,IJ.s) + c 8a (<j)I 8V {m,H2,nz) < 
4e-^ + 10- 7 e"i i£ ^A oo ( O -) + lO" 101 , 
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where the functions I pp , I ps , I sp , and I ss are evaluated at Z := 2(^2) if prj < "0 an d a t Z := max j£{i,2}{£w(/j 2 ),Ai' (MM2))}> 
> °o- 

Proof: We use a computer to calculate r\p(pi, .Z) for i G {1, 2, 3} and we prove that these quantities are bigger than 

> Z, what 



1. As z(a) < Z (see 11.27), rip(p,i,Z) > 1 for i £ {1,2,3} implies that |£| < r^^ 2 and that TV iAt3 



simplifies I ss (see equation (11.35)). We estimate the integrals as it is shown in the appendix using a computer, 



taking 5q = 1 if [i\mv > 10 and Sq = jk if p\mv < 10. We use the computer again to show that (6.12) is valid with 



(4e 2a 'i + 10 _3 e~i3 (M2 ?" ) A°(/i 2 )) instead of (4e 



10^ 3 e-l i;II ¥ ai A (CT)), (4e + 10- 7 e~I i±i ^A oo (^ 2 )) 



instead of (4e 



10 7 e 34 2 A°°(ct)), — Z instead of £ and CT(/ii) instead of ct(ct) (for T <E {pp,ps, sp, ss}). 



Finally by Remark 6.2 and the fact that ct(c) < ct(p-i), It(pi, P2, fJ-3,0 < It(pi, P2, P3, — Z), T e {pp,ps, sp, ss} 



(see (11.27)) , we obtain (6.12). 



□ 



6.1 The Incoming Electron Wave Packet. Continued 

THEOREM 6.6. Suppose that the magnets and energies are the ones of the experiments of Tonomura et ai. Then 
for every gaussian wave function with variance a € [— , ^] and every C, £ K with £ < —z(a) we have, 



\\ e - z ^ H W-(A)ip v - xe~ % ™ Ha <Pv\\ < e 



S;e{i, 1/2,0, -1/2,-1} °°cr l + 10 



-420 



(6.13) 



where the quantities A, t 00 are explicit numbers that depend only on the magnet and the energy that we take (see (4-32)) 



Proof: Equation (6.13) is a consequence of Theorem 3.1 Remark 4.1 and Lemma 6.1 



6.2 The Interacting Electron Wave Packet. Continued 

THEOREM 6.7. Suppose that the magnets and energies are the ones of the experiments of Tonomura et ai. Then 
for every gaussian wave function with variance a € [^p, ^] and every ( £ R with \(\ < z(a) we have, 



\e- t ^ H W-{A)(p v - xe^/o °° <r - A ( x + TV ) dT e - l ^i Ho ip v \\ < 



(6.14) 



4e 2^ 



E, 6{ i,i/ 2 ,o,-i/2,-i}(l + 10- 3 )A^ + 10 



-101 



10 



-420 



where the quantities A? are explicit numbers that depend only on the magnet and the energy that we take (see (4-32)) 



Proof: Let a € [— , y], then there are /zi, /z 2 and /Z3 such that fix, (12, P3 and a satisfies the hypothesis of Lemma 6.5 



We prove using a computer that they satisfy also the hypothesis of the Theorem 3.5 We obtain (6.14) from Theorem 
13.51 Remark 14. II and Lemmata|6.1l 16.51 



□ 
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6.3 Outgoing Electron Wave Packet and Scattering Operator. Continued 



THEOREM 6.8. Suppose that the magnets and energies are the ones of the experiments of Tonomura et ai. Then, 
for every gaussian wave function with variance a € [— , -y] and every C, £ K with C > z (o~) we have, 



We'^W-iA)^ - e l *xe~ 4 ™ ff Vv|| < 



7e -^ + e -I^ E« 6 u,i/2,o,-i/2,-il(l + 10- 7 )Af ^ + 1(T 101 + 1(T 420 , 



||5y> v - e**xyv|| < 



7e -^ + e-I 1 ^ Ei e{ i,i/ 2 ,o,-i/2,- 1} (l + 10- 7 )Af <r* + 1(T 101 + 1(T 420 , 



(6.15) 



(6.16) 



where the quantities Af° are explicit numbers that depend only on the magnet and the energy that we take (see (^.32)). 



Proof: Let a G Then, there are fJti,fi2 and /i3 such that /Lti,/i2,M3 an d cr satisfies the hypothesis of Lemma 



6.5 We prove using a computer that they satisfy also the hypothesis of the Theorem 3.8 We obtain (6.14) from 



Theorem 3.8 Remark 4.1 and Lemmata 6.1 6.5 To get equation (6.16) we remember that to obtain the error bound 



in Theorem 3.8 we used the error bound for the scattering operator of Theorem 3.7 Then, the error bound that we 



get for the outgoing wave function in Theorem 3.8 bounds the error bound for the scattering operator. 



7 Aharonov-Bohm Ansatz. Discontinuous Change of Gauge Formula 
from the Zero Vector Potential 

In this section we denote by A the vector potential constructed in Section |4j We take also the parameters, magnets 
and energies introduced in Section [5] 

7.1 Statement of the Aharonov-Bohm Ansatz 

Let Ax and A 2 be two differentiable magnetic potentials defined in K 3 \ K with curl zero and that have the same flux 
$. Suppose, furthermore, that 



\Ai{x)\ < C i+ , <z;(r) := rasx x&t3 ^ x ^ r {\Ai(x) ■ x\} € I/^O.oo). 
Choose any point xq € K \ K. We define 



\a 2 ,a 1 {x) ■= f {M - Ax), 

Jxq 



(7.1) 



(7.2) 



where the integral is over any curve in M 3 \ K that connects Xq with x. This integral does not depends on the curve 
because both potentials have curl zero, and both have the same flux <I>. If this last condition is not true we can not 
define A^ 2 a . Then, 
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A 2 =A 1 + WX A2 ,a 1 - (7.3) 

The solution to the Schrodinger equation with magnetic potential A 2 and initial condition given when the time is zero 
by the estate ip, is obtained in terms of the corresponding one for the magnetic potential A\, by the change of gauge 
formula, 

e -iiH(A 2 )^ = e i\ A2 , Ale - l iH(A 1 ) e ^\ A2 , Al ^^ 

The solution to the Schrodinger equation for the vector potential A x that behaves as 

e- l i Ho ^ (7.5) 



when the time goes to minus infinity is given by the formula (see equation 1.9 1, 



-HH^W-iA^ip. (7.6) 



In other words, ( 7.6 ) is the solution to the Schrodinger equation when the initial conditions are taken at time minus 



infinity by (7.5). Now we give the change of gauge formula for the Schrodinger equation with initial conditions taken 



at time minus infinity: 

e-W^W-iA^ip = e iA ^,A 1 W e -^ff(Ai) w -_( Ai ) e -iA A2 , Al , oo: (-p)^ ) ( 7J ) 

where Xa 2 ,Ai,oo( x ) : ~ um r-s-oo ^A 2 ,A 1 ( r x). (see equation (5.8) in [3]). 

Although the magnetic potential, A, constructed in Section [4] has curl equal zero, it has non zero flux. Therefore, 
there is no change of gauge between the vector potential zero and A. Suppose now that for every time the electron is 
practically localized in a region, T>, that has no holes (that is simply connected) or, in other words, in a region where 



Xa,o can be defined by equation ( 7.2 ) if we take curves that connects Xq with x lying on this region. On this region A 



is gauge equivalent to the vector potential zero and the change of gauge formulae (7.4) should follow approximately 
(although not exactly, because there is not a real change of gauge between A and the zero potential). The error will 
depend on how much of the electron lies in the complement of V. This is the Ansatz of Aharonov and Bohm [3J. Let 
us be more specific. In our case we take, 



V := (M 3 \ K)\S, (7.8) 



where 



S := {( Xl , x 2 , 0) e M J : yj x\ + x\ > f 2 }. (7.9) 
For two vector potentials A\ and A 2 whose curl is zero (and that do not necessarily have the same flux) we define 



the function given in (7.2) in the simply connected region T>: given x$ = (afyi, xq^ 2 , Xq^) € T> with xq^ < —h and x 
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in T> we define, 

Xa 2 , Ai (x) := / (A 2 - A x ), (7.10) 



where the integral is over any curve in T> connecting x with x. Note that for an electron to cross from the negative 
vertical axis to the positive one over X>, it has to go through the hole of the magnet. 

Then, we have that, 

A 2 (x)=A 1 (x)+WX A2tAl {x), xeV. (7.11) 

We extend \a 2 ,Ai to 1R 3 \ K by zero without changing notation, i.e., \a 2 ,a 1 {x) = 0, for x e S. Note that Xa 2 ,a 1 is 
discontinuous on S. 

The Ansatz of Aharonov and Bohm can be stated in the following way. 
DEFINITION 7.1. Aharonov- Bohm Ansatz with Initial Condition at Zero 

Let Ai be a magnetic potential defined in M 3 \ K such curl A\ — 0, and with flux not necessarily zero. Let ip the 
initial data at time zero of a solution to the Schrodinger equation that stays in T> for all times. Then, the change of 
gauge formula ([2], page 487), 

e -^ H ( A ^ip fts ip AB (x,t) := e i^ 1 ,oW e -«liJo e -iAA 1 ,oW^ {j 12) 

holds. 

Note that if the initial state at t = is taken as e~ tXA ^- a ^ tp the Aharonov-Bohm Ansatz is the multiplication of the 
free solution by the Dirac magnetic factor e lXA i- ^ [7]. 



Equation (7.12 1 is formulated when the initial conditions are taken at time zero. Now we reformulate it taking 
initial conditions when the time is minus infinity and for the high velocity state (p v . For the high- velocity state tp v 
and for big v, we have that, 



e 



For this statement see the proof of Theorem 5.7 of [3]- Formula (7.7) with W-(0) = I, and equation (7.13) suggest the 
following formulation of the Aharonov-Bohm Ansatz, with initial condition at time minus infinity and for high- velocity 
states. 

DEFINITION 7.2. Aharonov-Bohm Ansatz with Initial condition at — oo. General Potentials 

Let A\ be a magnetic potential defined in ]R 3 \!iT such curb4i = 0, and with flux not necessarily zero. Let ip v (Ai)(x, t), 

MAi){x,t) ^e-'^^WLCAi)^ 
be the solution to the Schrodinger equation that behaves as 

e-'i fl V v (7.14) 
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when the time goes to minus infinity. We suppose that tp v (Ai)(x, t) is approximately localized in T> for every time. 
Then, the following change of gauge formula follows, 



i/j v (Ai)(x,t) « e ttA I ,oW e -i|flo e -ttx 1 ,o, „(-v) v)v) 



(7.15) 



where o,oo(ac) = lim r -+ 00 \A 1 fl(rx). 



a 



Let us show that formula (7.15) can formally be derived from (7.12). We take -0 = e*^ A i'°^e l *Ai,o,<»( v >ip v m 



(7.12). Then, we have that e~ l s H ( Al V ~ e iAA i-°(*)e _i 'R- ffo e _<XA i' '°° (_ *Vv- For big velocities, the time evolution 



e l fi Ha ip v is localized near the classical position v£ [8j . Therefore, 



and thus, e 's 5 ^ 1 ^ behaves as (7.14) when the time goes to minus infinity. Then, 



and (7.15) follows 



For a general C 1 vector potential with curl equal zero and flux $, there is a real change of gauge (given by 



formula (7.2)) between this potential and the vector potential A with support in the convex hull of K constructed in 
Section |4j As the vector potentials A and A\ are gauge equivalent, they define the same physics and, therefore, we 
can always chose the vector potential A. For this potential, \a,o,co{~ v) = 0, and then, the Aharonov-Bohm Ansatz 
for initial conditions at minus infinity and the potential A is as follows. 

DEFINITION 7.3. Aharonov-Bohm Ansatz 

Let A be the magnetic potential constructed in Section Q). 
Let ip v (x,t) :— e~ % Ti H ( A } W— (A) ip v be the solution to the Schrddinger equation that behaves like 

VVo := e" i|ff Vv (7.16) 

when time goes to minus infinity. We suppose that ip v is approximately localized in T> for all times. Then, the following 
change of gauge formula holds, 

W « i>AB, v {x, t) := ^.'We-'^Vv (7.17) 

□ 

Observe that the Aharonov-Bohm Ansatz is the multiplication of the free solution by the Dirac magnetic factor e lXA ' ^ 
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Note that as we noticed before, the electron -when emitted, would follow the free evolution e~ l i Ha ip v under 
the assumption that we take a representation where the magnetic potential (A) vanishes at this time. If we take a 
representation given by a general vector potential (A\) with flux $, we should change the initial conditions at minus 
infinity by e'^i^-'-^e-'^"^ (notice that \ Au o,oo(-v) = *Ai,A,oc(-v) ). 



In the following sections we give a rigorous proof that (7.17) holds and we obtain error bounds for the difference 



between the exact solution and the Aharonov-Bohm Ansatz. We also provide a physical interpretation of the error 
bound and we relate it to the probability for the electron to be outside the region V. 

8 The Time Evolution of the electron Wave Packet. Final Estimates 

In this Section we use the same symbol, e~ l ^K H °, for the restriction of the free evolution to A and, moreover, we 
designate by |j • || the norm in L 2 (A). 

8.1 Incoming Electron Wave Packet. Final Estimates 



LEMMA 8.1. For every gaussian wave function, if, with variance a and for every £ € 
following estimate holds. 



10 



-420 



with £ < — z(o~), the 



(8.1) 



Proof: Let V-h be the set {(xi, x%, £3) G M : X3 < —h}. We have that, Xa,o{x) — and x( x ) = 1 for 1 £ 'D-h- Using 



polar coordinates we obtain (see (3.14), (11.3) and Remark 11.1) 

4 



lX ^.o e - l ^ H "ip v \\ 2 < 



r 3/2 



;3YD_„-vC)pKC) 



(8.2) 



Finally we notice that V2e 2 < 10 _4d4 for a > <tq. 



a 



Using Theorem 6.6 and Lemma 8.1 we prove that, 



THEOREM 8.2. Aharonov-Bohm Ansatz. Incoming Wave Packet 

Suppose that the magnets and energies are the ones of the experiments of Tonomura et al. . Then for every gaussian 
wave function with variance a € [ — , ^f] and every £ € R with £ < —z(a), the solution to the Schrddinger equation 



that behc 



7.16) when the time goes to minus infinity, e t ^i H W-(A)ip v , is given at the time t = ^ ( £ being th 



- C 



vertical coordinate) by the Aharonov-Bohm Ansatz, 

e iA A ,o(z) e -i-|ff 

up to an error bound of the form: 



<f v , 



(8.3) 
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\e^ i ^ H W-(A)(p v - e lX ^ a e- l ^ Ha ip v \\ < e-i 12 ^ 



(^{1,1/2,0,-1/2,-1}^°°^' + ^) + 10 



-419 



(8.4) 



where the quantities A i 00 are explicit numbers that depend only on the magnet and the energy that we take (see 



8.2 Interacting Electron Wave Packet. Final Estimates 

LEMMA 8.3. For every gaussian wave function, tp, with variance a € [A.b/mv, fi/2] and for every ( 6 1 with 
ICI < z {°~)> the following estimate holds. 



Hxe^io °° *-^(z+Tv)dT e -i4// 0(/?v - e iXA -°e' i ^i Ho ip v \\ < 2e-5 r iP( CT <C) 2 < 2.0031 e"^ + 



(8.5) 



2e -I^ +10" 456 . 



Proof: We denote by HM := {(x 1 ,x 2 ,x 3 ) G M 3 : y/x\ + x\ < n}. For x e HM, - f Q °° v ■ A(x + tv)g?t = A^ (x) 
and x(a;) — 1- Using polar coordinates we obtain (see ( |3.14| |11.3| )), 

4 



\ xe -*fa °°^ A ^+T^dT e -i^H 0(pv _ e iAA,o e -i4^o v , v ||2 < 



(R 3 \«M-VC)P(^:C) 



e" x < 4e- r i p( ^ c) . (8.6) 



The second inequality in (8.5) is proved in three cases: 



By (6.9 1, see also Sections [2] and 5.2 



(8.7) 



cr e [cr , 3.2 x 10~ 6 ]. For these values of a we have that tD(cr) = 2000~ 1/2 . We use |6^2| ), ( |ll.23| and the triangle 
inequality for the square-root term to obtain, 



z(a) < 1 (68h + V2000 x 34 



o~ 2 mv mv \ a 2 a 



Then, 



exp 



2(T 2 1 j_ 1 <■ 68/i , V2000x34 n 2 



exp 



2 ct2 +(^f(^ + V2000 x 34)2 



(8.9) 



The function f(a) = 1/ (a 2 + ^2 (^ + V2000 x 34) 2 ) restricted to the interval [cr , 10~ 7 ] has derivative 
equal to zero on the positive axis only at the unique point of intersection of the function cr 4 and the line 
j~fi(68h + V2000 x 34cr), see Sections [2] and 5~2 For the interval [10~ 7 , 3.2 x 10~ 6 ] the derivative of / is zero 
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over the positive axis in the unique solution of the equation a 4 = ^^2(68/1+ (V2000 x 34 + 680)cr), see Sections 
[2]and[!r2] Then, it follows that, 



exp 



2 cr 2 + ( S§^ + V2000 X 34) 2 



< 



(8.10) 



maxue{v , 10- 7 , 3.2xio- 6 } ex P 



2 ! ,2 + _l^(SM + y2000x34) 2 



Evaluating (8.10 1 using the experimental energies and magnets, we find that, 



e -ip( CT ,C) 2 < 10 - 



45* 



.11) 



• a€ [3.2 x 10- 6 ,^] 



Now we use that 



r 2 lP ^-0 2 < e -^ e -(rlp{cr,z(<r)) 2 -^) = 



exp 



2 (■ *Q) )2 
1 v a 2 mv t 

^ i + (4^) 2 



.12) 



By ( 11.23 ) g ^ is decreasing as a function of cr (see Sections |2l and 



5.2 



and notice that t™") _ i s decreasing 



on a) and then, we have that, 



V2 



exp 



2(T 2 1 4. f )2 



< \/2 



exp 



■ z(3.2xl0~ 6 ) x2 

■ (3.2xl0- 6 ) 2 mt;/ 



2(3.2 x 10 



-6^2 



, / z(3.2xl0-a) . 2 
1 V(3.2xl0- 6 ) 2 ml;' , 



< 1.4171. 



(8.13) 



REMARK 8.4. The term appearing in the middle inequality of equation (8.5) is two times the square root of the 



probability for the free particle to be outside the hole of the magnet (HM) when the electron is classically at the 
position (0, 0, £): 



E 3 \«M 



\{e- l ^ Ho ip v )(x)\ 2 dx = e ~ r * p{a ' c)2 



.14) 



Recall that HM is defined in the proof of Lemma 8.3 Equation ( 8.14 ) is a measure of the part of the electron that hits 
the magnet when the classical electron (the electron under classical mechanics rules) lies within a distance less than 



z(a) from the center of the magnet. By the second inequality in (8.5 ) we can see that the probability of the electron 

1 r2 

to be outside the hole of the magnet at time Q/v splits in two terms: one, e~ 2 , is due to the probability of the free 



electron to be outside the hole when C = (see formula (8.14)). This factor provides us an idea of the influence of the 
magnet over the electron given by the size of the wave packet (i.e., how much does the electron hits the magnet -see 
Section 



33 (mu)^ 

9.4), and the other, e _ 53 2 , is related with the spreading of the electron as time increases - see Section 



9.5 



This factor is important when a is small, because by Heisenberg uncertainly principle when the electron is localized in 
a small region, its momentum is not localized and therefore the electron spreads. Those two factors are essentially the 
causes of all the error bounds that we have in this paper. The error bounds are mainly produced by the probability of 
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the electron to hit the magnet when it is classically at the position (0, 0, (), with |£| < z(a). In Section[9]we provide 
an analysis of these terms and we give precise definitions of the size of the electron wave packet and of the opening 
angle, that is due to the spreading. 

Using Theorem |6 . 7| and Lemma |8.3| we prove, 
THEOREM 8.5. Aharonov-Bohm Ansatz. Interacting Electron Wave Packet 

Suppose that the magnets and energies are the ones of the experiments of Tonomura et al. . Then, for every gaussian 
wave function with variance a € ' iH an< ^ ever V C G ^ with |£| < z(a) the solution to the Schrddinger equation, 



1 vh H W~{A)ip v , that behaves as (7.16) when time goes to minus infinity is given at the time t = ^ ( £ being the 



i 



vertical coordinate) by the Aharonov-Bohm Ansatz, 



up to an error bound of the form: 



\\e-' l ^ H W-(A)(p v - e^.oe-^^Vvll < 

( 8 - 16 ) 

6.0031e"^ + e-l^^Eie/i i/ 2 0.-1/2 -n(l + 10- 3 )A<V + 2) + lO" 101 + 1(T 420 + 1(T 456 , 



where the quantities A? are explicit numbers that depend only on the magnet and the energy that we take (see (4-32)). 

8.3 Outgoing Electron Wave Packet and Scattering Operator. Final Estimates 

LEMMA 8.6. For every gaussian wave function, ip, with variance a and for every Q € K with ( > z{o~), the following 
estimate holds. 

Hxe^e-^Vv - e tXA -°e-^ H ^ v \\ < V^e''^^ 1 + 10" 420 . (8.17) 

Proof: Let T>h be the set {{xi, X2, £3) £ M 3 : x$ > h}, note that Xa,q(x) = 4> and x( x ) — 1 f° r x <= ^h- The proof 
follows in the same way as the proof of Lemma |8.1[ 



□ 

Theorem |6.8| and Lcmma|8.6| imply the following theorem. 



THEOREM 8.7. Aharonov-Bohm Ansatz. Outgoing Electron Wave Packet 

Suppose that the magnets and energies are the ones of the experiments of Tonomura et al. . Then for every gaussian 
wave function with variance a € [-^j , 9*] and every ( £ R with ( > z{a) the solution to the Schrddinger equation, 



e l vh H W~(A)tp v , that behaves as (7.16) when the time goes to minus infinity is given at the time t = - ((( being the 



vertical coordinate) by the Aharonov-Bohm Ansatz, 

e iXA '°e-^ Ho (p v , (8.18) 
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up to an error bound of the form: 



2 
1 1 



(8.19) 



7e -^ + e-^E i6{lil/2i0rl/2i _ 1} (l + l(T 7 )AfV + y/2) + 10- 101 + 2 x 10- 
and, furthermore, the scattering operator satisfies, 
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HS^v-e'Vvll < 7e~^ + e -l i£ ^( ^ (1 + l(T 7 )AfV + %/2) + lO" 101 + 2 x 1(T 420 , (8.20) 

1/2,0,-1/2,-1} 



where the quantities A°° are explicit numbers that depend only on the magnet and the energy that we take (see (4-32)). 



8.4 Uniform in Time Estimates for the Electron Wave Packet 



REMARK 8.8. The error bound of Theorem 18. 21 is smaller that the one of Theorem 18.51 and this last one is bounded 
by the error bound of Theorem |8.7| This is physically reasonable, because for an electron to be an interacting electron, 
it has to be first incoming electron and for an electron to be outgoing electron it has to be before an interacting electron, 
so the error should be accumulative. Let us prove this. That the error of Theorem |8.2| is smaller than the one of the 



Theorem 8.5 follows directly from the definitions (4.32). To prove that the error in Theorem 8.7 bounds the one of 



Theorem 8.5 we use again ( |4.32[ ) and that (remember that amy > 4.5), 



(1 + 10- 3 )^ + (2- y/2) = (1 + 10- 3 ) [^r^(l + 1.11 x 10- 6 )V2i|m^ + (2 - y/2) < 

(2 - 10- 3 )(150(1 - h)o-mv - 134.99)| ^— ^m 5 < (2 - 10~ 3 ) (a 1 / 2 mvn - a - V 2 i34.99/i)^M < ( 8 - 21 ) 



(2 - 10- 3 )(a 1 / 2 A^ + a-^ 2 AZu 2 ). 



□ 



This gives us the following theorem. 
THEOREM 8.9. Aharonov-Bohm Ansatz. Time-Uniform Estimates 

Suppose that the magnets and energies are the ones of the experiments of Tonomura et al. . Then, for every gaussian 
wave function with variance a € [rpj and every £ 6 K the solution to the Schrodinger equation, e~ l: ^ H W-(A)ip v , 



that behaves as 



by the Aharonov-Bohm Ansatz 



7.16) when the time goes to minus infinity is given at the time t = ^ ( £ being the vertical coordinate) 



.22) 



up to an error bound of the form: 



\\ e - l ^ H W-(A)ip v - e^.oe-'^Vvll < 



7e -^ + e-l i ^(E Je{1 , 1/2 ,o.-i/ 2 .-i } (l + 10- 7 )AfV + y/2) + lO" 101 + 2 x 10 



-420 



.23) 
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Moreover, the scattering operator satisfies 
HSVv-e^vll < 7e~^ 



33 (<rmu)- 
g 34 2 



^ (1 + 10" 7 )AfV + y/2) + 10~ 101 + 2 x 10" 

1/2, 0,-1/2,-1} 



120 



.24) 



XTie quantities A°° are explicit numbers that depend only on the magnet and the energy that we take (see (4-32)) 



By (4.29) and ( 5.2 1 m t (xi) > ^-ife), i € {1, • • • ,5}, and as a ll ' 1 mvr 1 > 



134.99 h 

a 1 ! 2 



□ 



(crmfri > 134.99 h, remember 



that amv > 4.5), we have that A J (a, v, m(xi)) > A J (er, v,m(x2)) (see (4.33)). We have also (see (4.32) and (4.33)) 



that A°°(cr, Vi, m(xi)) < A°°(cr, Vi,m(xi)) for i € {1,2,3} and j € {1,2} (notice that A(m)i > _4(m)4 and ^(771)3 > 



^4(771)5). So if we write A°°(wi, m(xi)) in (8.23 8.24) instead of Af° we obtain also error bounds, but now the 



coefficients A°° are fixed for all the magnets and velocities. Taking this into consideration we calculate the values of 
A («i, m(xi)) and we obtain the following theorem. 

THEOREM 8.10. Aharonov-Bohm Ansatz and Tonomura et al. Experiments 

Suppose that the magnets and energies are the ones of the experiments of Tonomura et al. . Then, for every gaussian 
wave function with variance a € [— , and every (G M, the solution to the Schrddinger equation, e~ t: ^ H W-(A)(p v , 



that behaves 



7.16) when the time goes to minus infinity is given at the time t — £ ( £ being the vertical coordinate) 



by the Aharonov-Bohm Ansatz, 



.25) 



10- 



2 x 10" 



up to an error bound of the form: 

\\e- l -^ H W-{A)ip v - e lXA - e- l ^ Ho tp v \\ < 

r 2 2 

7e - ^ +e-l if! ^(1.04 x 10 14 cr + 3.91 x 10 8 cr 1/2 - 1.41 x 10 3 - 1.14 x lCT 2 ^ 
Furthermore, the scattering operator satisfies, 

r 2 2 

||SVv - e i<& ^ v || < 7e~^ +e-I izl ^(1.04 x 10 14 cr + 3.91 x lOV 1 / 2 - 1.41 x 10 3 - 1.14 x lO^^H 
10~ 101 +2x HT 420 . 



.26) 



.27) 



□ 



We now bound the right hand side of 



.26)by7e ^ + ^(a, m«), where F(a, mv) := e ™ ( " (1.04 x 10 14 er 



3.91 x 10 8 cr 1 / 2 ) + 10 101 + 2xl0 420 . We notice that T is decreasing for mv fixed and amv > 4.5. We compute J~(15. 5 x 
10~ 10 , 1111)3) and show that this quantity is less than 10~ lo °, it follows that T{p, mv) < 10 -100 for a > 15.5 x 10~ 10 and 
the experimental velocities. Then T{a, mv) < e"I i£Z ^ (1.04x 10 14 (15.5x 10~ 10 )+3.91 x 10 8 (15.5x 10" 10 ) 1 / 2 ) + 10" 100 
= 177 x 10 3 e _ tf " 2 + 10~ 100 . We obtain the following theorem, that is our main result, and that is quoted as 
Theorem 1 1.1 1 in the introduction. 
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THEOREM 8.11. Aharonov-Bohm Ansatz and Tonomura et al. Experiments. Final Estimates 
Suppose that the magnets and energies are the ones of the experiments of Tonomura et al. . Then, for every gaussian 
wave function with variance a £ [— , fy] and every £ £ R, the solution to the Schrddinger equation, e~ r ^ H W-(A)(p v , 



that behc 



L mv ' 2 

_ c 



by the Aharonov-Bohm Ansatz. 



7.16) when the time goes to minus infinity is given at the time t = ~ (C, being the vertical coordinate) 

e lXA - o{x} e-^ H "<p Vl (8.28) 



up to an error bound of the form: 

We-^W-iA)^ - e iXA <°e- l ^ Ho tp v \\ < 

7e~^ + 177 x lO^-l 1 ^ + 10- 100 . 
Furthermore, the scattering operator satisfies, 

\\Sip v - e*Vll < 7e-^ + 177 x lO^-^^ 1 + 10" 100 . (8.30) 

REMARK 8.12. In the experiments of Tonomura et al. [26], they send an electron wave packet that partially hits 
the magnet. The part of the electron wave packet that hits the magnet does not go behind the magnet because we 
can see the black shadow of the magnet behind it. In other words, this part of the electron wave packet will be in the 
region {(x\, X2, £3) £ A : 13 < h}. We can bound, therefore, the probability of interaction of the electron with the 
magnet by the probability for the electron to not be behind the magnet for large time. We denote, as in the proof 



of Lemma 8.6 by T>h the set {(x\, xi, x%) g I 3 : 13 > h}. Actually T>h is the region behind the magnet. Then the 



probability of interaction of the electron with the magnet is bounded by, 

\\XA\v h e- lL - H W-{A)vA 2 (8.31) 

when the time goes to 00, where XA\v h is the characteristic function of the set A\T>h- We take as before ( = vt, then 
we have, 

\\XA\v h e~^ H W-(A)cp v r < (We-^W-iA)^ - e ^o e -i^H ^| + \\ XAK v h e-^ H °<p v \\) 2 . (8.32) 

We take Lu(cr) := j= and z{a) := z&i a -\ a (h(a)), see Sectional Using polar coordinates we obtain for £ > z(a) 

(see Section M (|3.14[), (|11.3|) and Remark [iTT]) , 



\XA\V h 



r 3/2 



e x dx < -e 



1 33 /„ m „A2 

— e _ 3i( amv > 



(M?\v h -vOp(°,0 



.33) 



Letting the time go to 00 in (8.31) and using Theorem 8.11 ( 8.32 1 and (8.33) we obtain that the probability of 
interaction of the electron with the magnet is bounded by, 



(7e - ^ + 177001e-i if ^ 



10- 



.34) 
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9 Physical Interpretation of the Error Bounds 



We analyze the error bounds given in equation (8.29 8.30). The error bounds appearing in the whole paper are pro- 



duced by the same factors. Equations (8.29 8.30) provide uniform in time error bounds that apply to all experimental 



magnets and energies. The behaviour of the error bound is the same for the three energies and the two magnets, so 
there is no loss of generality if we select a magnet and an energy in our analysis. We will use the biggest energy (E\) 
and the second magnet {K%) to provide numbers and graphics. So, for now on we take the magnet K-i and the energy 
Ex. 



The main factors that produce the error bound in equation (8.29 8.30) are the terms 



1. Size of electron wave packet factor. 



2. Opening angle of the electron wave packet factor. 



e 2» 



33 (Tmri- 

e 34 2 



(9.1) 



(9.2) 



When the variance a is close to the radius of the magnet, (9.1) is close to 1 and (9.2) is extremely small, because 



in this case amv is big. Then, when the electron is big compared to the inner radius, (9.1) is the important term, 



which justifies our name. When the variance is very small -such that amv is close to 1- the factor (9.2) is close to one 



and (9.1) is extremely small ( ^ is big), and then, the important factor is (9.2). But when the variance in position 
(a) is small, by Heisenberg uncertainly principle the variance in momentum is big, and then, the component of the 
momentum transversal to the axis of the magnet is large. In consequence, the opening angle of the electron wave 



packet is large, and the electron spreads fast as it propagates. This justifies the name of (9.2) 



By the previous discussion, we divide the analysis of the error bounds in (8.29 8.30) in three sections: big sigma 



(cr close to the inner radius of the magnet), small sigma (amv close to 1) and intermediate sigma (sigma neither big, 
nor small). 

9.1 Big Sigma, <re[|,|] 



Remember that n = f i — e and that e is defined in Section 5.2 Here f\ — 1.75 x 10 cm (see Section |5.l|. Then, 



in terms of absolute values, big sigma ranges over the interval [7.7955 x 10 -6 , 8.7500 x 10~ 5 ]. In Figure [T] we show 
the graphic of the en 
representative values. 



the graphic of the error bound in (8.29) as a function of ~, for big sigma, and in the table below we give some 
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Error Bound as a Function of 
Sigma Over r% for Big Sigma. 


Sigma Over t\ 


Error Bound 


.34305 


io- 1 


.27626 


io- 2 


.23764 


io- ;i 


.21170 


10" 4 


.19274 


IO" 5 


.17811 


10~ 6 


.16637 


io- 7 


.15668 


IO"* 


.14851 


10" y 


.14150 


10 -1(J 



9.2 Intermediate Sigma, o G [6.7591 x 10 6 r 1; ||], or a G [- 



23 1546781 



l mti 1 mv 



Remember that mv — 1.9842 x 10 10 (see Section 5.1 ). Therefore, in terms of absolute values, intermediate sigma ranges 
over the interval [1.1592 x 10~ 9 , 7.7955 x IO -6 ]. For these values of sigma, — > 22 and amv > 23, and therefore, the 



error bound in (8.291 is less than 10 



-99 



For intermediate sigma the probability of interaction of the electron with the magnet is less than 10 (see Remark 



.12) 



r 4.5 23 1 
L mv ' mv J 



9.3 Small Sigma, o e [1.3224 x 10"%, 6.7591 x lCTS], or a e [; 
In terms of absolute values we have that a G [2.2679 x 10" 10 , 1.1592 x 10~ 9 ]. In Figure |] we show the graphic of the 



error bound in (8.291 as a function of ■jf-, for small sigma, and in the table below we give some representative values. 



Error Bound as a Function of 
Sigma Over r\ for Small Sigma. 


Sigma Over n 


Error Bound 


1.6001 xlO" 6 


10" 1 


1.7234 xl0~ e 


10" 2 


1.8384 xlO" 6 


10" 3 


1.9467 xl0- B 


10" 4 


2.0492 xl0~ 6 


IO" 5 


2.1469 xl0" e 


10" B 


2.2403 xlO" 6 


io- 7 


2.3299e xlO" 6 


10" b 


2.4162 xl0~ 6 


io- y 


2.4996 xl0- B 


io- 10 



9.4 The Radius of the Electron Wave Packet 



As before, we denote by T-LM the cylinder {(x, y, z) € M 3 : \Jx 2 + y 2 < ri}. MM is basically the hole of the magnet. 
The factor e^S 5 is practically the square root of the probability for the free particle at time zero to be outside the 
hole of the magnet: 
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e 20 



Xm\umv 



e 2" 



(9.3) 



This factor represents the part of the electron wave packet that hits the magnet or goes outside (the square root 
appears because our estimations are in norm and not in probability). It is natural to have this factor in the error 
bound because we are only modeling the particles that go trough the hole. This factor is significant only when the 
variance is close to the inner radius of the magnet. As the proximity of the electron to the magnet increases the error 



in equations (8.29 8.301, it is important to define intuitively what is the meaning of this closeness or, in other words, 
what is the size of the electron wave packet. We agree that the free electron is actually localized in configuration space 
in a ball centered in the classical position vi and with radius chosen in such a way that the probability of finding the 
electron on this ball is 99%. We measure the radius of the wave packet at the time t = - when the free particle is in 
the center of the magnet - and denote it by R(a). Then, we have: 



R := R(a) = 2.382 a. 



The error due to the part of the electron that hits the magnet (9.3) is practically zero (smaller than 10 99 ) when 
R < .10827-1 (i? < 1.8556 x 10~ 5 ). In Figure |] we show the error bound of equation (|8.29[) as a function of the radius 
of the wave packet over n for big sigma, a £ 9-] (.1082 r\ < R < .5102 n). 

Even when the size of the wave packet is comparable to the inner radius of the magnet we have error bounds 
extremely small. We give some data to show this behavior: 



Error Bound as a Function of the Radius 
of the Wave Packet Over r\ for Big Sigma. 


Radius of the Wave Packet over r\ 


Error Bound 


.81716 


10" 1 


.65806 


10" 2 


.56606 


10~ 3 


.50427 


io- 4 


.45911 


10" 5 


.42425 


10" e 


.39629 


IO" 7 


.37322 


10"" 


.35376 


io- y 


.33703 


io- lu 



9.5 The Opening Angle of the Electron Wave Packet 



Although it is impossible to define an opening angle of the electron, because it is everywhere, we agree to say that the 
free electron (in momentum representation) is actually in a ball, Bp(Mv) with center the classical momentum (Mv) 
and radius P such that there is a 99% probability for the electron to have its momentum within this ball. We define 
the opening angle, oj(a), in the obvious way (see Figure |4j), 
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. ,w(<t), P 2.382 

sm( ) := = . 

2 Mv amv 

When sigma is big, the opening angle is very small and when sigma is small, the opening angle is big, this is 
nothing more than Heisenberg uncertainty principle. 



The factor e 34 ( 2 ' of the error bound 



.29 



30 1 has the following interpretation in terms of the opening angle: 



e 34 2 

This factor is practically zero (smaller than 10 -100 ) when oj < 11.8 degrees (<r > 1.1592 x 10 -9 or amv > 23), and 
then, it begins to increase as lj increases ( a decreases). In Figure [5] we show the error bound in equation ( 8.29| as a 
function of the opening angle for small sigma, a £ [1.3224 x 10~ 6 ?'i, 6.7591 x 10~ 6 ri], and in the table below we give 
some representative values. 



Error Bound as a Function of 


the Opening Angle for Small Sigma. 


Opening Angle (degrees) 


Error Bound 


51.8407 


10- 1 


47.8885 


10~ 2 


44.7231 


10~ 3 


42.1135 


10- 4 


39.9137 


10~ 5 


38.0265 


10" e 


36.3842 


10- 7 


34.9380 


10- 8 


33.6517 


10- 9 


32.4979 


10 -io 



10 Conclusions 



In Theorems |8.2[ |8.5[ |8.7[ |8.9[ |8.10| and |8.11| we found the time evolution of the electron up to an error bound that 
we provide explicitly. The approximate wave function of the electron that we give is the one given by the Aharonov- 
Bohm Ansatz. It coincides also with the part of the electron wave packet that goes through the hole of the magnet 



in Tonomura et al. experiments |26j. As we noticed before (see Section 7.1) the Aharonov-Bohm Ansatz is valid if 
the evolution of the exact wave packet is localized at every time in a simply connected region, with no holes, (for 



example in (7.8)). The main factors that produce the error bounds are the size of the wave packet (see (9.1)) and 



the opening angle (see (9.2)). These factors can be understood also in terms of the part of the wave packet that 



hits the magnet when the electron crosses the hole of the magnet (see Remark 8.4) and, therefore, they are related 
with the part of the electron not localized in a simple connected region (see (7.8)) at every time. In Section [9] we 
analyzed the error bounds and we have shown that our estimates for the time evolution are valid for a rather big 
interval that starts when the opening angle is close to 55 degrees ( a ~ 1.3224 x 10~ 6 ri ) and ends when the size of 
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the wave packet is close to the inner radius of the magnet (close to ri). We have shown also that the error bounds 
decrease very fast -exponentially- as the variance gets away from the extremes of the interval. For intermediate sigma 



(cr € [6.7591 x 10 °r%, ||]), the time evolution given by the Aharonov-Bohm Ansatz (8.28) differs from the exact one 
only by a number less than 10~" in norm. As it is shown in Remark 8.12 and Section 9.2 for intermediate sigma, 



the probability that the electron wave packet interacts with the magnet is smaller than 10 199 and so, there are no 
fields in the trajectory of the electron. Nevertheless, the solution is the one given by the Aharonov-Bohm Ansatz 



(8.28) and it is affected by the vector potential A by a phase factor e lXA -° . This phase factor is the one that appears 
in Tonomura et al. experiments |26j . Although in the experiments of Tonomura et al. |26] there is no interaction with 
the magnetic field, there is an interaction with the impenetrable magnet. Tonomura et al. |26j argued that it is not 
necessary to consider the part of the electron wave packet that hits the magnet -they used a rather big one- because 
the shadow of the magnet was clearly seen in the hologram. Our results show that it would be quite interesting to 
perform an experiment with a medium size electron wave packet with an intermediate sigma. One could use, as well, 
a bigger magnet. Our results show that quantum mechanics predicts in this case the interference patterns observed 
by Tonomura et al. |26j with extraordinary precision. 

In the Aharonov-Bohm Ansatz the electron is not accelerated, it propagates following the free evolution, with the 
wave function multiplied by a phase. As we prove that the Aharonov-Bohm Ansatz approximates the exact solution 
with an error bound uniform in time that can be smaller that 10 in norm, we rigorously prove that quantum 
mechanics predicts that no force acts on the electron, in agreement with the experimental results of Caprez et al. [6]. 

Summing up, the experiments of Tonomura et al. |17| I25[ 126] give a strong evidence of the existence of the 
interference fringes predicted by Franz [9^ and by Aharonov and Bohm [2]. The experiment of Caprez et al. [6 verifies 
that the interference fringes are not due to a force acting on the electron, and the results of this paper rigorously 
prove that quantum mechanics theoretically predicts the observations of these experiments in a extremely precise way. 
This gives a firm experimental and theoretical basis to the existence of the Aharonov-Bohm effect 2j, namely, that 
magnetic fields act at a distance on charged particles, even if they are identically zero in the space accessible to the 
particles, and that this action at a distance is carried by the circulation of the magnetic potential, what gives magnetic 
potentials a real physical significance. 

11 Appendix A. Estimates for the Free Evolution of gaussian States 

In this appendix we prove estimates for the solutions to the boosted free Schrodinger equation, 

d 

i~^-<p(x,z) = Hi<p(x,z), ip(x,0) = (p(x), (11.1) 
az 



where the boosted free Hamiltonian H\ is defined in (3.11) 



Recall that under the change of variable t :— z/v, the solutions of ( |11.1| are solutions of the boosted free Schrodinger 
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equation with Hamiltonian e~ mv ' x H e lmv x . Classically, a particle that starts at the origin with velocity v = (0, 0, v), 
will be located at time t at the position (0, 0, z). At the high- velocity limit, the quantum evolution follows the classical 
one and the parameter z can be taken as the position in the z— direction of the particle. We consider the case where 
the initial state is gaussian, 

^-P^e^/- 2 , (11.2) 
with variance a. The solution to (JTOJ is given by, 

3/2 i 

e~ lzHl (fi — e -' lzmv / 2 - - e -(x-zir) 2 /2(a 2 +iz/mv) _ (11.3) 

7r 3 / 4 (a 2 + iz/mv) 3 / 2 

we will often use the following simple result. 

REMARK 11.1. Suppose that C 3 < C 2 < C x < 0. Then, 



Ci rC 2 -Ci roo 

e~ z dz < e~ c i / e- z dz < e" Cl / e~ z dz. (11.4) 

C 3 Jc 3 -Cx Jo 



° 2 z 2 e- z2 dz < -%-° 2 + %~ c * + \e~ c * Cl e~''dz < e^-^ 2 - + \ I e^' dz). (11.5) 

c 3 2 2 2 Jcs-d ^ 2 J 



Proof: 



° 2 e~ z2 dz < e-° 2 f 2 e -(*"-C?) < e " c ? f^J-^'dz, 



e 1 e 

C*3 J C3 



where we used that, z — Cf > (z — C\) . This proves 1. Furthermore, 2 follows from 1 and the following equation. 



2ze~ z dz = --e~ z 2 + / -e~ z dz. 



c 3 2 2 7 C3 2 



LEMMA 11.2. Xef f be a bounded complex valued function with support contained in D. Then, for z > h and 
d > h — z, 



\f(x)e- lzH ^\\ < e-'W-.*) 2 ^ 
V2 



v2 



(11.6) 



(11.7) 



/(x) e -""v|| < e"*-^ 2 / 2 ^ r 2 p(z) 3/2 . (11. 
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Proof: We use the function p{z) denned in (2.9), 



-izH lr J\ 2 <r 11/11 



7T 3/2 J(D-vz)p(z) 



dx < 



11/11' 



2 r (h-z)p(z) 



,1/2 



(-fo-z)p(z) 



dp e^ 2 ( 



1 - e- r ^ 



< 



,1/2 



oo e -e« n „( CT , 2 ) 2 



-2/ip(z) 



(11.9) 



where in the last inequality we used (11.4). Equation (11.6) follows from (11.9). Equation (11.7) is obtained similarly. 



Equation (11.8) follows from (11.9) and the estimate, 

rO 

-27ip(z) 



dze-* 2 (l ~ e^pM 2 ) < 2hr 2 p(zf 



LEMMA 11.3. Let f be a bounded complex valued function with support contained in D. Then, for Z > h, s > 0, 



/~ \\f(x)e- i ^\\ < 



%/2 e 



/ 2 (max(Z,s) - Z)+ 



ll/IU -8j„„(j,max(Z,s)) 2 /2 



(11.10) 



2hr 2 {amv) 3/2 / t 



max(Z,s) ( ( r 4 m 2 t ) 2 +C 2 ) 3 / 4 



Proof: We prove the lemma writing the integral in the left hand side of (11.10) as follows 



oo />max(Z,s) 
izHi 



||/(x)e- l *"V|| = 



||/(x)e-^V|| + / H/^e-^VH 

J max(Z,s) 

and using (11.6) in the first integral, (11.8) in the second, and the fact that 9i nv (a, z) 2 is increasing in z for z > h. 



□ 



LEMMA 11.4. Let g 



p3 , ,p3 



&e bounded and with support contained in D and let z > h. Then, 



\g(x) ' P' 



-izH 1 II < Ijgjloo -e mv (a,z) 2 /2 



,(o-,z) 30T 



1/2 



(11.11) 



5 (x)-pe- l ^V < 



llffll 



00 



rV4 f 



'^' 2 [A(amv) 2 + 2] 1/2 Vhr 2 p{zf' 2 . 



(11.12) 



Proof: 



\g(x)-pe-" H ^\r < 



7r 3 / 2 o- 2 



da; < 



J^[T(a, z) + 6(a, z)] (l - e -^ 2 ) , (11.13) 



/(D-vz)p(z) 

where and T are defined in Section [2] Equation ( jll.ll | follows from (11.13) applying the last inequality in (11.4) 
to T(cr, z) and the last inequality in ( |11.5| to 0(er, z). Furthermore, using the middle inequality in (11.5) we obtain 
that, 



(11.14) 



2hp(z) 
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Note that, 



/ e" z2 dz < 2hp(z). 

J-2hn(z) 



Writing z — h — z + h — 2h in (11.14) we obtain that 



Q(a,z) < ^(z + h)e- 0i ™^ 2 Azhp{zf < e~ 0i ^ {a ' z) \hp{z){amv) 2 . 



Moreover, applying the middle inequality in (11.4) to T(cr, z) we prove that, 

T(<r,z) < e- e -"^ z)2 2hp(z). 



(11.15) 
(11.16) 

(11.17) 
(11.18) 



Equation (11.12) follows from (11.13 11.17 11.18) 



LEMMA 11.5. Let g : R 3 — > C 3 be bounded and with support contained in D. Then, for any Z > h with 6i nv (a, Z) > 
I, s > 0, 



/~ \\g(x) -pe-^vll < |f|j| e -^^^ 2 /2 ( l*W^)l 1/2 + (max(Z, S )-Z) + 

W^ e -e I „„( CT ,ma X (z, s) ) 2 /2 ( 2ffmB + ^2) Vh r 2 {amvf' 2 .C^^m V + C 2 )~ 3/4 - 



(11.19) 



Proof: We split the integral in the left hand side of (11.19) as follows 



poo (>max(Z,s) poo 

/ \\g(x)-pe- izH ^\\= / || 5 (aO-pe-" ff V||+ / || 5 (x) • pe^V) 

JZ J2 Jmax(Z,s) 



and using (11.11) in the first integral, (11.12) in the second, and the fact that the functions e x / 2 */x,e x / 2 x 1 / 4 are 
decreasing for x > 1 (notice also that Oi nv (cr, z) 2 is increasing in z for z > /i). 

REMARK 11.6. Suppose that z, ( £ R + , s and 6 are real numbers such that z > £, s > z — 2£, & > 0. Then, 
1. In any interval I :— [<7i,<72] suc/i that's/a G J, — 6i nv (o~, z, s, £) > y/l/2, 



T(a, z, s, C)e' b P^ 2 < maxpT^, z, s, C)e- b "^ 2 , T(a 2l z, s, CK^ 2 ^ 



£. In any interval I := [<ti, 02] suc/i f/ia£ Vcr G /, — 6i nv (cr, z, s, £) > y/3/2, 



9(a, z, s, C ) e ~M^) 2 < m ax[e(ai, z, s, C)e~ MCTl ' z)2 , e(<7 2 , z, s, Q e - bp ^ 3 ' t) ' 



Proof: We give the proof of 1. The proof of 2 is similar. We have that 



(11.20) 



(11.21) 
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T(a, z, s, C) e"M^) 2 = \ (g Cw) ((^f - **) \e~^^ (C + s)p(a, z) - 

(11.22) 

e-^-^ 2 P^ 2 (z-C)p(a,z) -2bp{a,z) 2 T(a,z,s,C) ■ 



As the function e x x is decreasing for x > 1/V%, the term in the square brackets is (11.221 is negative. Then, the 



left-hand side of (11.22) is different from zero for a S I if \J zjmv ^ / and otherwise, it is negative for a < \J z/r 



and it is positive for a > \f zjmv. This proves 1. 



□ 



Remember that z^^ih) is defined in Section [2] It is given by, 

h(amv) 2 <rmv (_ 9 9 (amv) 2 W^ 2 

zc.Ah) = - ( \ 2 ~_ 2 +77 , 2 -- 2U /2 "~V + h 2 ? \ 2 -_ 2 -1 • (11-23 

REMARK 11.7. Suppose that a 2 <G <a 1 . Then, 

zq,*{() < max(zu ttT1 (C),^, CT2 (C))- (11.24) 



Proof: Note that as a function of cr, p(cr, z) is increasing for er < \J zjmv and that it is decreasing for a > \J z/r 
Suppose that za,, ai (C,) < z Q.o-(C)- Then, a < \J zjmv, because if a > \J z/mv, 



dT 1 = -t 



■inv 



(o-l, Zq >(Ti (C), ZQ,ai(0>0 < -9inv((T, 2<J,<7i(0> ^,<ti(0>0> 



since, —9i nv {o~, z, z, () = (z — (,)p(o~, z) and as —9i nv is increasing in z > 0, this implies that ^,<t(C) < z <i,o-i (0- Then, 
02 < cr < -y/ z/mv, and it follows that, 

But as also, 

and —9inv is increasing in z > 0, we have that £aj )CT (0 < z^. (T2 (0- 

LEMMA 11.8. Let /Xj, z € {1,2,3} belong to R + . Suppose that the following conditions are satisfied, 

1. Either pi < o-q, * € {1, 2, 3}, or pi > era, i <E {1, 2, 3}. 
2- pi < P3,i € {1, 2}, or p { > p 3 , i € {1, 2}. 
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We define ^ raax := max(/i b /i 2 ), ^ min := min(/i 1 ,^ 2 ) and take v = ^ max , if ^ > fi 3 ,i £ {1,2} and v = ^ min , if 
Hi < M3,« € {1,2}. We denote by Z := z(^i max ), if m < er ; and Z := max ie(1 ,2}{^( Mmax ), w (%max))}, */ Mi > CT o- 



We suppose that Z > Vs.^WMniax))- Let f 



C be a complex valued function and we take 



fa,z{x) ■= f(x, er, z). Suppose the support of f a ,z is contained in K — [0, (z(c) — ( — z)]v for some ( 6 M, every a € R+ 
and every z € R z + £ < z(<r). Then, for every gaussian wave function ip with variance a € [/i m in , Mmax] , 



|/^(x) e -^e-^V|| < ^A»(A*i,Ma,Ai3,C) 



(11.25) 



where, 



/ ps (/i 1 ,/i 2)M 3,C) ■= ^^^(M^^eO.. iM2} e-^"^'^.^ 3 ( fc ^» )» 2 +7TV4 max {-C,0} 



(11.26) 



Proof: It follows from equation (11.231 that z(<r) < z&ro-Vo-CMMmax))- If Mi ^ °o then cl>(cr) = <l>(/i max ) for <r G 

[Mmin 5 Mm 



It follows from Remark 



11.7 



that z(o) < max ie { li2 }{^(^ max ), Mi (/i(Mma X ))} = 2. If Mi < o"o then from 
formula ( |11.23| and the definition of Q{a) we have that ^(^^(MMmax)) < z ii>( Mmax ), AImax (^(Mmax)) = z(Mmax). We 
conclude that 

z(cr) < Z, (11.27) 

and then, 



z((T)-C 



dz \\f (r<z (x)e 



—izBi—iCH-i 



Z-C 



<P|| < / dz ||/ (T , z (a;)e" " "■ ! « 
o 



(11.28) 



As in (11.9) we prove that 



\\LA^y- izHi e- tCHi 

Then, 

J Q Z ~ C \\faAx)e- izHl e- iCHl v\\ ^ < ft 



ll/l 



2 

OO 



V|| < ^lyf T(tr, z + C, Z, fc(/w)) e ~ r ^'^ . 

7T ' 



c(/° T(a, z, Z, /i(Mmax)) 1/2 e-^^^ 2 , 0)+ 



Jo 



^.^C'Omax)) 



(11.29) 



T(a, z, Z, /i(/w)) 1/2 e-^^) 2 + f s 



< 



I. 



7rV4 max (~C, 0)e- Irt'-O 2 + vr 1 ^ ^ {h{^))e-^ p ^^ ^ 
T(a,z,Z,/,( Mmax )) 1 /2 e -^P(^) 2 



(11.30) 



V2,*,M3 ( fe (Mmax)) 



where we used that, T < y 7 ^. If z > w ^ 3 (MMmax))j it follows from Remark 11.7 that z > CT (MMmax)) f° r every 



cr belonging to the interval limited by v and /i 3 . We complete the proof of the lemma using (11.20) in the integral in 



the right-hand side of (11.30), and for the other two terms we argue as in the proof of (11.20) 



□ 
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Using the proof of the preceding lemma, we prove the following, 

LEMMA 11.9. Suppose that the hypothesis of the Lemma \ll.8\ are fulfilled and furthermore, assume that the support 

of fa,z is contained in K for every a £ K+ and every Then, for every £ 6 K with |£| < z{a), 

C \\UA^- izHl e- iCHl 4 < M^2J pp ( / x 1 , / x 2 , / x 3 ), (11.31) 

where, 

Ipp{^l,^2, M3) := Ips(Hl, M2, M3, 0). (11.32) 

and 



^ \\U, z (x)e- izHl <P\\ ^ ^fe/ PP (m,M2,M3). (11-33) 

7T ' 



LEMMA 11.10. Let /j, max , fj, m i„,iy and Z be as in Lemma 11.8 We suppose furthermore that Z ^ z rj {h([i max )) 

V 3 '^M3 



and r\p{[ii, „ ^ 3 (h(p max ))) > 1 /or i G {1, 2}. Lei : K 3 x R + xM->C 3 &e a complex vector valued function and 

we take g a ,z{x) '■= g{x, a, z). Suppose that the support of g a ^ z is contained in K — [0, z(a) — C ~ z ] v f or some 

all a £ K+ and for all z with z + £ < z{a). Then, for every gaussian wave function ip with variance a £ \p m in> l^max], 



r (<7 C ||g^(^- P e-^ e ^V|| «fa< JIfJ^/ ss (/ii,/i2,M3,C) (11-34) 
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where, 



1/4 -^piUi,* fj (hfamax))) 2 

7 ss (Mi,M2,M3,C) := ^^^ i ^ 3 WMma X ))max M . G{w!M3} (e )- 



7T 1 / 4 



max 



{-C,0}max Mie{pi!M2} (e-^^^ 2 ) + 



tt 1 / 4 max(-C, 0) 



71-1/4 VWa^^^^elw.M^iKMi,^ 

iax ^e{Mi^2}( r iP(^,C)e~^ p(AllX)2 ), if |CI <V,^ 
e" 1/2 , */ |CI > r„ llWl 

^ 1/4 ^^(M^))m^ €{B , M (e^*-^- 3 (%«))) ! ) + 
tt 1 ^ max (_ C , o) max Mie{Ml)M2} (e- ! ?' , ^'« 2 )+ 

E we{ ^ 3} //_ (M/w)) e( W ,r ) Z > ft( AW )) 1/2 e-^^ a dr+ 



Proo/: By (11.27) 



o Jo 



Estimating as in the proof of (11.13) we prove that 



We have that, 



11^,, (a:) • pe-^e-K^H 2 < [&(*, z + (,Z, /i(/w)) + 

T(a, z + C,Z, /i(/w))(l + r\p{a, z + () 2 )] e -M^+C) 2 . 



where, 



Jl := m aa; (/ c (e(a, r, Z, ftf^)) e^*" 
(Mc»»)) 



-»^ 1/2 dr,0) + 



Jo' 



e -r lP ( CT ,r) ^ 



^6(ct, t, Z, /i(/x max )) « 

1 /2 

(e(a, r, Z, fc(/w)) e" r2 "^ T ) 2 J dr, 
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1 /2 

h ■= max{J° (X((7, t, Z, Maw)) rjp(tr, T )2 e -^,r)^ dT; ) 
+ Jo^""" 3 (M " max)) (T(a, r, Z, fc^)) r?p(*, r) 2 e -^(^) 2 ) V * dr, 

max (^y2. l -, M3 (' 1 (A , max)), min(r„, M3 ,Z)) ,,\ 1/2 

(T(<7,T,Z ! / l (M roax ))r-?p(a,T) 2 e-'-^^T) \ ^ 



/ 4 := 



T(a, t, Z, /z(/w)) r 2 p(a, r ) 2 e -^(<^) ) dr, 

"»l*V5,i/ 1 „l»K l «»«i 1 ™l r »,l«i''« 



1 /2 

7 6 := max(/° (t(<t, r, Z, ft(/w)) e^^) 2 ) dr,0) + 
£W*"--)> (T(a,r,Z,h(^))e-^fy /2 dr , 



z 



Since T < y^r and < t/it/2 we have that, 



/*::= / (T( ( T,r,Z,/.( Almax ))e- r2 ^ 2 ) ' ,/r 



1/2 



/i +/ 6 < ^ 1/4 max(-C,0)e-^K^O + ^^^(,w)) e 
■fTrV* Lax(-C, 0) e - 4p(-0 2 + ^ e~ * '<"'*>/Ws <M/w))) 2 A 

By Remark |ll.6| 



V 3'"^3 



(e( M i,r, Z,/.( Aimax ))e- r2 ^^ 



J 2 < J! / ((-)(//,. r.Z./M/w)W '"'") Jr. 



1/2 



/X < e{^^ 3 } ,/z V5,v, ff3 (fc(^max)) 

Moreover, since xe _:E / 2 is increasing for < x < 1 and decreasing for x > 1, 



h < T 1/4 max(-CO) I r M^C)e ; P( °'° > */ ICI < W a I 
I e-V>, i/|C|>W 2 J 



1/2 

dr. 



By Remark 11.6 if r > ^^^(^(Mmax)), 



T((T,t, Z,/i(/Xmax)) < max T(fl i ,T,Z,h(fi max )). 
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By (11.501 and as le x takes its maximum at x = 1, 



f Z 1/2 



(11.51) 



Note that if r\p{pi,T) 2 > l,pi £ {v,pz} then, r 2 p(a, r) 2 > l,Vcr between and ^3 . Hence, as in the proof of 



Remark 11.6 we prove that, 



< 



max MiS{ ^ 3} rfp(/i i ,r) 2 T(/x i ,r,Z,/i(M max ))e r ' p(Al, ' T)2 Xn Mie{ „, M3} K/^Cmlt)^!}^)) 



(11.52) 



and then, 



7 4 < 2 



max ( z ^2,- #3 (' 1 (Cm.x)),min(r„ >(13 ,Z)) 



1 /2 

(r 2 / 9(/i,,r) 2 T(^,r,Z,/ l ( / i max )) e -'-i^ I *' T ) 2 ) dr. (11.53) 



To obtain equation (11.34), we use (11.38 11.46 11.49 11.51 11.53) and we argue as in the proofs of Remark 11.6 



to estimate equations (11.46) and (11.49) 



□ 



Using the proof of the preceding lemma we prove the following, 



LEMMA 11.11. Suppose that the hypothesis of the Lemma \ll.lO\ are fulfilled, assume furthermore, that the support 
of g a ^ z is contained in K , for all a € R+ and for all z. Then, for every ( el with |£| < z(o~) and every gaussian 
wave function <p with variance a € \p m ini f-max], 



^ \\gaAx)-Pc- tzHl e- lCHl 4 dz<m^2I sp (p 1 , f i 2 , f i 3 ), 



(11.54) 



wht 



J"sp(Ml:M2,M3) := 7^(^1,^2,^3,0). (11.55) 

And 

\" || 5CT , 2 (x)-pe-^V|| dz< ^§^I sv {iix,^%,^). (11.56) 

LEMMA 11.12. Let f : R 3 — > C be bounded and with support contained in D. Then, for Z > h, and £ such that 

Z C \\f(x +(Z~(z + C))v)e-^ e -< ff V|| dz < (Z - £)J!/lb e -i fl '-«(*.^ a . (11.57) 

v 2 



!><S 



Proof: Estimating as in the proof of(11.9) we prove that, 

\\f(x + (Z-(z + C))v)e-^|| 2 < Mfr *~* *** Io 2P(a,Z+ ° e-^^dr 



(11.58) 



< 



-6 irlv (a,Z) 2 



where we used (11.4) 



LEMMA 11.13. Let g : R 3 — > C 3 be bounded and with support contained in D, suppose that 8i nv (a, Z) 2 > |. Then, 
for Z > h, and £ such that C < Z, we have that, 



\g(x + (Z -(z + C))^)-pe- tzHl e- lCHl ip\\ dz < (Z - 



\\g\\ 



,(cr, Z) Sa/tt 



1/2 



(11.59) 



Proof: The lemma is proven estimating as in the proof of (11.11 1 using Remark [TTTTJ 



12 Appendix B. Upper Bounds for the Integrals 



In this appendix we prove upper bounds for the integrals appearing in the terms I ps , I pp , I ss and I sp (see (11.26), 

( TQ2| ), ( JiQB] , ( fHT55l ». 

Suppose that Z > s > C, > 0. Designate := {0, 1, y/2, \/3, • • •}. We denote, 

(12.1) 



{Zi,Z 2 ,"' ,Z K } ■= V^ovNn [-6 inv (a,s,s,C),-d inv (a,Z, Z,Q], 
where Z\ < Zi < ■ ■ ■ < Zk- As — 0i nv (a, r, r, £) is increasing as a fuction of r we have that, 

LEMMA 12.1. Suppose that Z > s > r > 0, and let f : R ->■ K satis/?/ /(r) > r - 2(. T/ien, 
Jf dTT(<j,r,f(r),C) 1/2 < 



x / 4 r -l^ 
e 2 C 



7r 1/a r 



(a,s,s,Cy 



C*V,„(0 - ») + T,fJi e~W (z zr i (0 - z z -iJ0) + e- L * z « (z - % - v (C))] , 

1 ' J j + 1 > j ' k ' 



r V4 -^p(z iJO) 2 _1, 
V2 & 



e 2 



(Z Z 7\A0~S) + 



K-l -^2^ z z 



tM(*„-i JO) 2 



(%- „(0-%rv(C)) +e~^^) e-* z *(Z-z z -i (0)], 



2 r 



(12.2) 



(12.3) 



(12.4) 



A'-l -TPl J Z 



JO) 2 _i 



1 



(Z i+1 + V^/2) 1/2 (%- , CT (C) - %r\ CT (0) + 



(12.5) 



e -* P (z) e 4z t ^. w((7) Zj Z; C) + V5F/2)V2 (z _ ^ (c)) j 
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If moreover, r±p(Z) > 1, 



,.2 



// drT{a,Tj{T)X) 1/2 r lP {T)e-^^) < ^ [r lP (z z -, a (()) 
e -**T^ e -¥^M? [z A0 _ s) + rip[z „(0)e- ip{Z * 

1 ' J J + 1 ■ 



. „ (C)) 2 1 7 2 



(12.6) 



Proof: We split the integral in the left-hand side of (12.31 as follows, 



// drT((T,r,/(r),0 1/2 =j;* rl " (0 drT ( CT ' T > /( T )> C) 1/2 + T,f=i L^TjO drT{a^ f{rU) l/2 
f z _ x io drT(a,Tj(r),0 1/2 , 



(12.7) 



and we apply (11.4 1. This proves (12.31. (12.41 is proved in a similar way. Equation (12.51 is proven in the same way, 



but using (11.5). Finally, we prove (12.6) as above, using (11.4) and observing that the function xe x I 2 is decreasing 
for x > 1. 
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Error Bound as a Function of Sigma over n 
for Big Sigma. 
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Figure 1: Error bound as a function of a over r\ for big sigma 
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Error Bound as a Function of Sigma over n 
for Small Sigma. 
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Figure 2: Error bound as a function of — x 10 6 for small sigma 
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Error Bound as a Function of the Radius 
of the Wave Packet over r, 
for Big Sigma. 
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Figure 3: Error bound as a function of the radius of the wave packet over r\ for big sigma 
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Figure 4: Opening angle 
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Error Bound as a Function of the Opening Angle 
for Small Sigma. 




Figure 5: Error bound as a function of the opening angle for small sigma 
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